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Foreword 


This  Final  Scientific  Report  covers  the  work  performed  under 
Contract  AF  61(052)-916  during  the  period  from  1  January  1968 
to  31  March  1968  which  terminates  this  program.  It  is  pub¬ 
lished  for  information  only  and  does  not  necessarily  represent 
the  recommendations,  conclusions  or  approval  of  the  Air  Force. 

This  contract  with  the  "Laboratorium  flir  Werkzeugmasohinen  und 
Betriebslehre  der  Technischen  Hochsohuie  Aachen"  was  initiated 
under  Manufacturing  Methods  Project  "Chatter  Behavior  of  Heavy 
Machine  Tools".  It  has  been  accomplished  under  the  technical 
direction  of  Mr.  Floyd  L.  Whitney  of  the  Advanced  Fabrication 
Techniques  Branch,  MATF,  Manufacturing  Technology  Division,  Air 
Force  Materials  Laboratory,  Wright-Patterson  Air  Force  Base  OHIO. 

■»  i 

Closely  related  efforts  in  the  field  of  vibration  analysis  are 
covered  under  Contracts  AF  33  (61B)-2661  at  the  Cincinnati  Mil¬ 
ling  and  Grinding  Machines,  Inc.,  Cincinnati,  Ohio;  AF  33 
(6l5)-2664  at  the  Department  of  Meohaniaal  Engineering  Univer¬ 
sity  of  Cincinnati}  AF  61  <O52)-001  at  the  "Institut  ftlr  Werk¬ 
zeugmasohinen  und  Betriebswiaaensoh&f ten  der  Teohnisohen  Hoch¬ 
sohuie  Miinohen"  and  AF  61  (052)-966  at  the  "Laboratorium  fUr 
Werkzeugmasohinen  und  Betriebslehre  der  Teohnisohen  Hoohsohule 
Aachen" . 

The  principal  investigator  of  this  project  Is  Professor  Dr.-Ing. 
Dr.  h.c,  D.  Bo.  H,  Opitz,  Chief  engineer  is  Dr.-Ing.  A,  Mussen- 
brock.  The  work  on  the  project  was  performed  by  the  following 
research  group; 

Research  engineers  on  chatter  analysis  and  Investigations  of 
the  rusting  prooess  are  Dipl.-Ing.  H.  Fiubeth,  Dipl,~Ing.  F.  Ber- 
nardi  and  Dipl.-Ing.  K.  Beckenbauer. 

Research  engineers  on  structure  analysis  are  Dr.-Ing.  W.  Dapper 
and  Dlpl.-wirtaoh.-Zng.  II.  Groth. 
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Abstract 


This  project  has  been  engaged  in  a  program  of  research  which 
was  directed  at  the  problem  of  determining  and  improving  the 
dynamic  stability  of  machine  tool-metal  cutting  systems.  Such 
a  study  will  provide  machine  tool  manufacturers  with  the  neces¬ 
sary  analytical  methods,  tests  and  specification  techniques  to 
assure  a  chatter  free  cutting  process  for  large  machine  tools. 

Based  upon  the  theory  of  regenerative  chatter  the  system  ma- 
ahine-tool  -  cutting  process  is  reduced  to  a  nonintermeshed 
closed  control  loop,  of  which  the  stability  behavior  is  represen¬ 
tative  for  the  ohatterfree  cutting  performance.  For  economically 
carrying  out  the  stability  analysis  digital  computer  programs 
have  been  developed,  Results  of  practical  chatter  investiga¬ 
tions  are  finally  detailed. 

The  discussion  of  different  possibilities  to  provide  data  for 
chatter  spedif ioations  showed,  that  further  practical  investi¬ 
gations  must  be  carried  out,  before  one  of  the  disoussed  proce¬ 
dures  can  be  finally  stated  to  be  fully  suitable  and  before 
maohine  tool  chatter  specifications  can  be  published  respec¬ 
tively. 

Methods  and  digital  computer  programs  have  been  developed  and 
tested  for  the  precalculation  of  the  static  and  dynamic  pro¬ 
perties  of  maohine  tool  systems  already  in  the  status  of  design. 

A  short  general  view  over  this  research  work  is  presented. 

For  calculating  the  deformation  of  columns  with  square  and 
rectangular  cross-section  equations  have  been  derived.  They 
include  the  influence  of  shear  force  and  a  cross-sectional 
distortion  under  torsional  loads  can  be  considered  as  well. 

In  model  tests  these  equations  have  been  proved.  Furthermore 
the  influence  of  different  types  of  ribbing  on  the  static  and 
dynamic  behavior  of  columns  has  been  determined. 
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Introduction 


In  machining  processes  on  cutting  machine  tools  often  intense 
vibrations  occur,  which  can  shorten  the  tool  life,  increase 
the  wear  of  the  machine  and  cause  insufficient  accuracy  and 
surface  quality  of  the  workpieces.  Because  of  these  vibrations 
-  in  general  called  chatter  vibrations  -  in  many  oases  the  per¬ 
formance  of  the  machining  operation  aan  even  be  impossible 
without  changing  the  given  cutting  conditions.  In  the  plant 
a  change  of  the  cutting  conditions  is  not  necessarily  but  actu¬ 
ally  in  most  of  the  practical  cases  connected  with  a  decrease 
of  cutting  performance,  which  moans  a  forfeiture  of  profit 
for  the  machine  tool  user.  This  aspect  becomes  specially  im¬ 
portant  for  machines  which  require  high  amounts  of  investment, 
as  for  instance  large  and  heavy  machine  tools  and  NC  machines. 

The  above  listed  reasons  caused  the  United  States  Air  Foroe 
Materials  Laboratory,  Research  and  Technology  Division,  to 
sponsor  a  research  program  for  the  investigation  of  the  dynamic 
and  chatter  behavior  of  machine  tools  in  general. 

This  Final  Scientific  Report  covers  the  work  performed  under  Con¬ 
tract  AF  61  (052 )-916  at  the  Laboratorium  ftlr  Werkzeugmaschinen 
und  Betriebslehre  der  Technisohen  Hochsohule  Aachen,  The  re¬ 
search  work  was  initiated  on  1  January  1966  and  terminated  on 
31  March  1968.  This  contraot  is  the  continuation  of  the  research 
work  conducted  at  the  same  institute  under  Contraot  AF  61(052)- 
713  during  the  period  from  1  June  1963  to  30  November  1904, 

Under  the  preceding  contract,  which  for  the  first  time  dealt 
with  systematical  investigations  into  the  static  and  dynamic 
properties  and  the  chatter  behavior  of  heavy  machine  tools,  fun¬ 
damental  research  work  on  the  field  of  the  development  of 
measuring  equipment  and  test  procedures  as  well  as  initial 
theoretical  Investigations  concerning  the  chatter  behavior  in 
face  milling  operations  have  been  carried  out.  Additionally 

-  1  - 


n-« 

initial  stops  nave  noon  done  to  calculate  the  static  anti  dyna¬ 
mic  properties  of  machine  tool  elements  by  means  of  electroni¬ 
cal  computers. 

The  aim  of  this  new  contract,  in  which  the  important  results 
of  tho  preceding  contract  have  been  used ,  was  to  improve  the 
dynamic  and  chatter  behavior  of  cutting  machine  tools  in 
general  and  therewith  to  get  higher  profit  from  the  machine 
tool  capacity.  The  total  research  work  la  directed  at  impro¬ 
ving  the  behavior  of  actual  machines  in  the  shop  and  of  ma¬ 
chine  tools  being  In  the  status  of  design, 

The  object  of  the  first  section  "Chatter  Analysis"  was  to  inve¬ 
stigate  the  theoretical  relationships  in  maahine  tool  chatter 
and  to  get  knowledge  about  the  effect  of  the  machining  para¬ 
meters  and  structural  characteristics  upon  the  ehatter  behavior. 
Extensive  theoretical  and  experimental  investigations  have  been 
discussed  in  the  Quarterly  Technical  Reports  (see  Ref,  U.-lB), 
The  most  important  results  and  especially  the  practical  appli¬ 
cation  of  the  developed  procedures  are  detailed  in  the  first 
part  of  this  report.  In  this  part  a  mathematical  model  being 
generally  valid,  is  established,  by  means  of  which  the  stabi¬ 
lity  of  milling  and  turning  operations  can  be  analysed.  In 
this  model  certain  simplifications  of  the  exact  Mathematical 
relations  are  included,  which  however  have  negligible  influence 
upon  the  aoouracy  of  the  results  for  the  practically  interesting 
conditions.  These  simplifications  have  only  bejn  introduced 
in  order  to  make  the  chatter  investigation  to  be  carried  out 
easily  and  economically. 

One  essential  purpose  of  the  theoretical  chatter  turn  lysis 
Is  to  provide  tho  neaessary  data  for  chosing  the  machining  con¬ 
ditions  in  such  a  manner,  that  a  high  ohatter  froo  autting 
performance  can  be  reached.  In  order  to  take  into  uoaount 
and  to  vary  the  practically  interesting  parameters  easily  in 
the  theoretical  calculation  different  digital  computer  pro- 
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grama  have  been  developed,  by  means  of  whioh  the  stability  analy¬ 
sis  oan  bo  carried  out,  In  many  application  oases  a  very  good 
correlation  between  theoretioal  and  experimental  borderlines 
of  stability  has  been  obtained. 

An  additional  purpose  of  the  research  work  in  the  section 
Chatter  Analysis  was  to  study  the  basis  for  establishing  dyna¬ 
mic  specifications  for  maohine  tools,  Dynamic  (chatter)  speci¬ 
fications!  Which  are  demanded  for  sihoe  about  10  years  chiefly 
by  the  machine  tool  users,  shall  complete  the  geometrical  Spe¬ 
cifications  after  Sohlesinger,  Dynamic  specifications  Shall 
enable  maohine  tool  Users  to  , judge  and  to  compare  the  machines 
reliably  with  regard  to  Uioir  olmtterfree  cutting  porfrrmanoo, 

investigations  on  milling  machines  and  lathes  of  various  types 
and  manufacturers  gave  the  possibility  to  study  the  conditions 
for  the  establishment  of  general  dynamic  specifications,  From 
these  extensive  investigations  on  principle  different  possibi¬ 
lities  aan  be  derived  for  testing  ahd  specifying  maohine  tool 
chatter  behavior,  these  procedures  will  be  uisoussed  in  the 
second  section  "Chatter  Rpec if ioatiohs".  However,  considering 
the  numerous  aspects  concerning  chatter  specif  icmt  ions ,  whioh 
were  unknown  until  now  ahd  therefore  wore  not  taken  into 
account  when  initiating  this  program,  it  can  be  stated  here 
beforehand,  that  it  seems  Very  premature  to  establish  and 
to  publish  chattel*  spec  if  ioations  at  this  time, 

Another  part  of  this  contract  dealt  with  struotur-nnalysiai, 

From  the  above  mentioned  point  of  view  and  In  a  ',ot  of  other 
oases  it  would  be  very  helpful  for  surely  being  able  to  pre¬ 
calculate  the  static  and  dynamic  behavior  of  maohine  tool  parts 
and  systems  of  these  elements  whioh  consequently  give  a  ma¬ 
chine  tool , 

Within  this  part  cf  the  contract  Home  methods  wore  developed 
to  calculate  the  later  behavior  already  in  the  statue  of 
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design  (nee  Kef. [2  4  18j).  Starting  from  a  flux  of  foro©  analysis 
over  some  total  machine  constructions  the  Bteps  for  developing 
the  methods,  programing  and  testing  them  was  fixed. 

During  this  work  problems  occured  with  respect  to  an  econo¬ 
mical  preparation  of  input  data  as  for  example  cross  section 
values.  Por  thle  purpose  additionally  a  method©  was  found 
and  will  be  fully  described . 

The  methods  carried  out  here  may  be  called  solving  a  lot  of  so 
far  not  or  unsufficiently  solved  problems  by  aid  of  a  computer 
during  the  design  prooess,  The  usefulness  and  accuracy  of  the 
suggested  ways  are  shown  within  the  several  Quarterly  Techni¬ 
cal  Reports  comparing  measured  and  computed  results  about 
representative  structures, 

However,  it  has  to  be  pointed  out  already,  in  the  foreword 
that  Still  problems  Ooaur  When  applying  the  methods  to  arbitra¬ 
ry  chosen  systems  with  respeot  to  the  boundary  oondltiohs  in 
form  of  so  far  unknown  static  and  dynamic  behavior  of  the 
connection  points  between  elements  and  to  the  foundation. 

Resides  the  computation  of  the  static  ahd  dynamic  behavior  of 
machines  or  structural  elements  it  is  still  necessary 
to  test  actual  structures  themselves  or  models ,  In  order  to 
determine  the  Influence  of  internal  ribs  and  wall  ribs  on  the 
behavior  of  oolUmn  models  with  different  types  at  ribbing 
have  been  tented.  These  tests  point  out  some  interesting 
aspects  for  designing  columns.  In  addition  to  the  touts  an 
analysis  of  the  deformation  of  columns  loaded  by  a  single  force 
or  a  couple  of  forces  is  given  especially  with  respeot  to  the 
deformation  due  to  shear  force  and  aross-Heotlonal  distortion 
caused  by  torsional  loads, 
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1 .  Chatter  Analysis 

Vibrations  on  machine  tools  ooeuring  relatively  between  work- 
piece  and  tool  can  be  reduced  to  different  exciting  mechanisms 
which  however  in  a  certain  case  can  only  be  distinguished  hard¬ 
ly.  According  to  the  fundamental  research  work  into  the  vibra¬ 
tion  behavior  of  machine  tools  conducted  especially  by  Tobias 
and  Tlusty-Polacek,  the  observed  vibration  phenomena  can  be 
classified  into  two  main  groups: 

Forced  vibrations 

Self  excited  vibrations. 

Forced  vibrations  are  often  generated  by  external  forces,  which 
are  transferred  by  the  machine  foundation,  Unbalance,  errors 
and  defects  in  gear  drives  and  bearings  have  unfavourable  effect 
chiefly  in  these  cases,  where  their  frequency  correlates  to  a 
hatural  frequency  of  the  machine,  In  these  oases  resonances  are 
excited  with  often  very  great  vibration  amplitudes  at  the  cut¬ 
ting  point,  in  the  shop  the  disturbance  sources  in  general  can 
be  localized  and  removed  easyly. 

In  this  connection  the  periodia  force#  in  interrupted  cutting 
processes  and  especially  In  face  milling  by  the  tooth  mesh  im¬ 
pacts  as  well  as  by  the  varyihg  number  of  teeth  being  in  contact 
with  the  workpiece  are  of  high  importance,  Forced  vibrations,  as 
they  can  be  observed  in  face  milling)  are  characterized  by  the 
function  of  the  vibration  amplitude  at  the  cutting  point  versus 
the  deapth  of  cut,  which  is  given  in  Figure  1,1  a.  it  can  be 
seen,  that  until  a  certain  deapth  of  out  the  relationship 
between  vibration  amplitude  and  deapth  of  out  is  almost  linear. 
An  it  was  stated  in  practical  cutting  tests,  the  vibration 
amplitude!*  keep  almost  constant  with  Increasing  deapth  of  cut 
above  this  point.  In  the  plant  resonances  excited  by  periodic 
forces  in  general  can  be  avoided  by  suitably  choosing  the  cut¬ 
ting  conditions,  In  many  oases  of  milling  a  slight  change  of 
the  rotational  speed  is  already  successful. 
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However,  in  face  milling  operations  it  is  sometimes  difficult  to 
distinguish  forced  vibrations  from  self  exoited  vibrations 
immediately.  For  self  exoited  vibrations  the  function  vibration 
amplitude  versus  deapth  of  cut,  plotted  in  Figure  1.1  b  is 
characteristic .  It  can  be  seen,  that  the  vibration  amplitude 
increases  suddenly,  when  exceeding  a  certain  deapth  of  cut,  the 
socalled  critical  deapth  of  cut.  An  important  characteristic 
of  self  exoited  vibrations  is,  that  they  can  arise  due  to  a 
single  disturbance  and  then  continue  and  increase  respectively 
without  any  other  external  force  acting.  The  energy  being 
necessary  for  the  vibration,  is  here  produced  by  the  cutting 
process  itself.  In  thr  following  the  term  "Chatter”  is  exclusive¬ 
ly  reserved  for  vibrations,  which  occur  due  to  self-exa itation. 

Many  theories  have  been  developed  In  the  last  years  in  order  to  ex¬ 
plain  the  observed  self-oxcited  vibrations  by  different  physi¬ 
cal  effects. 

However,  as  it  was  found  out  by  Fishwick  Gurney,  Tobias 
[4,6)21]  and  Tlusty-Polaoek  [l]  ,  in  turning  processes  the  so- 
called  regenerative  effect  is  mainly  the  cause  for  chatter, 

Research  work  performed  by  P4ters-Vanherok  [20J  as  well  as  the 
invetitigations  performed  under  Contract  AF  01  (082)-7l3  [lo] 
showed,  that  the  regenerative  effeot  is  valid  in  milling  opera¬ 
tions  too, 

In  a  stationary  cutting  proaess  a  static  cutting  force  is  ac¬ 
ting,  which  onuses  a  static  deformation  of  the  elements  in  the 
l'lux  of  foroe,  If  this  static  force  is  superimposed  by  a  dyna¬ 
mic  force  -  for  example  an  initial  disturbance  due  to  material 
inhomogeneity  -  an  additional  relative  displacement  is  exoited 
between  workpiece  and  tool ,  which  produces  a  chip  thickness 
variation  and  a  wavy  workplace  surface.  The  amount  of  tho  dis¬ 
placement  and  thus  the  deapth  of  the  waves  is  dependent  on  the 
dynamio  compliance  of  the  machine.  If  the  tool  cuts  this  wavy 
surface  again,  a  dynamio  cutting  force  is  generated,  which 
acts  upon  the  machine  too.  This  mechanism  nun  bo  represented 
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by  the  block  diagram  in  Figure  1.2.  The  stability  of  this 
closed  control  loop  and  there  with  the  chatter  behavior  depends 
on  the  properties  of  the  machine  on  one  hand  and  on  different 
parameters  given  by  the  cutting  process  on  the  other  hand. 

For  the  stability  analysis  the  milling  process  was  at  first 
replaced  by  a  model  similar  to  that  used  already  for  turning 
processes  by  Toblaa ,  Tlusfcy  and  Pfeters.  Hereby  it  is  assumed, 
that  the  chip  thickness  variation  takes  place  in  the  middle 
of  the  arc  of  contact  and  that  the  cutting  forces  acting  at  the 
cutter  teeth  can  be  taken  as  one  resultant  force  acting  in  the 
middle  of  the  arc  of  contact  too.  With  this  assumptions  the 
machine  dynamics  can  be  measured  as  directed  response  locus, 
i.  e.  excitation  in  the  direction  of  the  resultant  force  and 
measurement  in  the  direction  of  the  maximum  chip  thiokness 
variation.  The  stability  analysis  on  the  basis  of  the  directed 
response  loaus  is  relatively  easy  to  carry  out,  however,  the 
directed  response  loaus  is  valid  for  limited  cutting  conditions 
only.  Another  decisive  disadvantage  Is,  that  the  approximation 
of  the  milling  process  by  this  so-called  single-tooth  model 
yields  the  more  incorrect  results  the  greater  the  arc  of  con¬ 
tact  Is, 

For  this  reason  the  multi-tooth  model  has  been  developed 
[8, 9,10, 11-14]  ,  whereby  the  orientation  of  the  chip  thiukuess 
-  and  cutting  force  variation  is  not  taken  into  account  al¬ 
ready  during  the  measurement  of  the  response  loci,  but  by 
additional  direction  coefficients  in  the  following  stability 
calculation.  By  this  the  dynamic  measurement  can  be  conduc¬ 
ted  completely  independent  from  a  given  machining  configuration. 
Usually  the  excitation  and  the  measurement,  aro  performed  in 
the  three  directions  of  a  rectangular  oar thou  inn  coordinate 
system,  of  which  the  axes  correlate  with  the  feed  directions 
of  the  machine,  By  introducing  the  direction  coefficients  now 
the  angular-position  of  eaoh  tooth  and  the  amount  and  the 
direction  of  the  cutting  forces  acting  at  each  tooth  can  be 


-  7  - 


noiiwH Herofl  _  ThU2 

proooss"  can  ho 
nans  coefficient 


f  Itti  n«in 


vitviu^vM  AM 


4.  U  ~  kS 
MIC  IJi.VUft 


II-..-*  *  4  .. 


U  U  I,  W  XII 


represented  by  two  single  blocks  "chip  thick- 
"  and  "direction  coefficient"  as  shown  in 


y; 


Figure  i,:t. 


Basud  upon  the*  theory  of  remunerative  chatter  iri  the  following 
at  first  the  system  equations  will  bo  formulated  and  the  stabi¬ 
lity  analysis  will  bo  explained,  Hornby  the  heavy  duty  machining 
of  Fe-matttrials  with  carbide  tools  shall  mainly  bo  considered. 

1  , 1  Control  Loop  Represent at  ion  of  t ho  Dynamic  Machining 
Process 


As  already  shown,  the  relationships  being  valid  in  regenerative 
chatter,  oan  be  represented  by  a  olosutl  control  loop.  The 
applicability  of  the  stability  criteria,  well  known  from  the 
oontrol  engineering,  depends  on  the  question,  in  what  mathe¬ 
matical  form  the  blocks  oan  be  described,  Therefore  it  seems 
reasonable  to  disauss  at  first  the  different  blocks  df  the 
control  loop. 

1,1.1  Discussion  of  the  Blocks  of  the  Control  Loop 

For  the  description  of  the  machine  tool  properties  with  regard 
to  a  chatter  investigation  it  is  sufficient  to  have  knowledge 
about  the  machine  dynamics  at  the  cutting  point  only.  Supposing, 
that  in  a  dynamic  machining  operation  dynamic  foraos  are  ac¬ 
ting  relatively  between  workpiece  and  tool,  which  can  bo  re¬ 
presented  by  their  components  !'x ,  Fy,  ¥z  in  throe  directions 
of  a  cartesian  co-ordinate  system,  the  dynnmia  behavior  at 
the  cutting  point  oan  be  described  b  nine  response  loot 
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whore  x,  y,  »  are  the 
between  worltpieco  and 

components  of  the 

t  ool  , 

relative  displacement 

Tho  response  loci  are  usually  measured  at  the  machine  itself , 
beoauso  up  to  this  time  no  sufficiently  exact  results  oftn  be 
obtained  using  digital  computation  methods.  For  this  purpose 
in  most  of  the  cases  exciter  teats  are  carried  out  at  the  ma¬ 
chine  a  tan ding  still,  whereby  static  and  dynamic  forces  are 

applied  between  workpiece  and  tool  by  means  of  special  exciters, 
The  relative  displacements  caused  by  the  dynamic!  forces  are 
rointed  to  these  forces  and  plotted  in  the*  forth  of  response  loci . 
If  an  automatic  transfer  function  analyv.er  ih  connection  with 
an  analbg-digitul-converter  and  a  tape  or  card  punch  is  used  for 
this  purpose,  the  data  are  available  for  the  stability  analysis 
on  the  digital  computer  immediately  after  the  dynamic  measurement. 

Under  Contract  Ay  01  (Ofl!2)-UU(i  on  the  basis  Of  t ho  random  noise 
theory  a  measuring  procedure  has  been  developed,  by  ntoaas  of 
which  the  dynamic  measurement  can  be  performed  at  the  running 
machine  and  oven  under  actual  machining  conditions  [ll»]  .  Good 
simulation  of  the  machining  conditions  is  possible  in  many 
cttMUK  especially  where  the  spindle  is  the  critical  struc¬ 
tural  element  -  using  the  contactless  uloatro-pmgnot io  exciter , 
developed  under  this  eon  tract;  (sou  [l?]  ). 

From  tho  properties  of  the  machine  the  data  chnraotori/.iug 
the  cutting  process  block  are  completely  independent.  Thu 
function  of  tho  Hooallod  direction  coefficients  cun  bo  ox-. 
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plained  looking  at  the  schematic  drawing  of  a  fade  milling  ope¬ 
ration  given  in  Figure,,  1^4 .  The  reiat,  lonahipB  derived  here  for 
milling  operations  are  also  valid  for  turning  operations,  how» 
ever  the  equations  describing  tho  turning  process  are  in  general 
more  simple, 

Assuming  at  firnt  a  dynamic  system  with  only  one  degree  of 
freedom  in  a-dircction  the  relation  between  a  relative  displace- 
ment  x  and  the  cutting  force  variation  at  the  i-th  edge 

of  the  cutter,  caused  by  this  displacement  OArt  be  represented  by 
the  diagram  in  'Figure  l,ft,  tho  cutting  force  variation  can  be 
formulated  as 

»  u(e"  -4)  h  -sin*  -Sinfi  -ccsJL  slnfftrfl)  (i>a) 

with  kQ  “  chip  thiokncdd  coefficient 

T  "  rotational.  titte  for  one  tooth  pitdh 

or 

Pi#  (k)  *  it’  '  W*  '^i'hk 

with  the  direction  coefficient 

*  d/'n  tf  •  sin  f;  '  cojjL  •  slnfft  +/3)  (1*4) 

Hereby  the  influence  of  the  inclination  angle  A  shall  bo  neglec¬ 
ted,  because  this  angle  is  practically  only  small, 

In  the  general  case  of  a  system  with  three  degrees  of  freedom 
tho  total  chip  thickness  variation  Uj  at  one  edge  results  from 
three  parts  u^x),  u^y).  u^Oi).  These  parts  are  caused  by 
relative  movements  in  x-,  y-.  and  m-dlrootion,  which  are  pro¬ 
duced  by  the  components  Flx,  F’iy,  of  a  dynamic  cutting  force, 
The  relationship!!  between  the  force  component*  and  the  compo¬ 
nents  of  tho  chip  thickness  variation  can  bo  derived  in  ana¬ 
logue  manner  as  the  equations  (1.2),  (l.H),  (l.*D.  Thus  the 
dirention  coefficient  are  obtained t 
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It  can  ho  seen,  that  .he  direction  coefficients  depend 

on  the 

out'tinii  force  direction,  tfiven  by  the  angles  and  /?  ,  on  the 
side  cutting  odi;e  angle  5«  and  on  the  angular  position^  of 
the  tooth.  Ah  the  angle  varies  continuously  against  the 
arc  of  contact,  tho  direction  aoef f iolwiitw  are  periodic  time- 
variant  functional  as  shown  in  Kl^um  1.(1.  In  milling  opera¬ 
tions  several  tooth  are  normally  in  oontnot  with  the  workpleow  | 
thus  the  (I iron t ion  coefl'ia lento  actually  are  the  sum  of  the 
values  being  valid  for  all  tooth t 
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Foi*m  (l.fi)  n.nri  (l.fl)  yields  the  equation  system  for  the  resul¬ 
tant  forces 

Fh  *  Ae  {Rxtt  UH  +  Rxy  Uy  f  RxzUZ^ 

Fy  kc  (  Ry«  ux+ fyy  Uz‘  (,J) 

Fz  •  kc  (  R?x  <vH  +  R*y  Uy  +  Hz*"*) 

It  one  oonslders,  that  a  chip  thickness  variation  can  aotually 
occur  at  those  teeth  only,  which  are  instantaneously  in  contact 
with  the  workpiece,  the  direction  coefficients  are  only  existing 
in  the  range  between  the  entrance  angle  Yg  and  the  exit  angle  ^ 
that  means  the  validity  of  the  equations  (1.5)  is  limited 


following  the  condition 

Rix*  ‘  1  ‘  Ritz  +  o 

f or 

r£  *  r 

>,  L  Y? 

1  "  'A  (l.S) 

Ri w*  •  •  *  Ri “0 

for 

Ta  *  Y 

’i  *  Te 

Because  of  the  Limitation  (1.8)  the  functions  of  the  direction 
coefficients  are  discontinuous  at  the  angles  of  tooth  entrance 
and  tooth  exit.  An  example  for  these  discontinuous  functions, 
calculated  for  u  given  milling  process, is  shown  in  Figure.  1.7, 

If  the  stability  analysis  shall  be  performed  applying  one  of  the 
stability  criteria  well  known  from  the  control  engineering 
technique,  it  is  very  time-oonsuming  and  there  with  not  econo¬ 
mical  to  take  into  account  tho  time-dependence  of  the  direc¬ 
tion  coefficients,  Therefore  extensive  fundamental  investiga¬ 
tions  have  boen  carried  out  [ 11, 12, 13]  ,  in  order  to  determine 
tho  uonditions,  under  which  the  timevariant  coefficients  aan  be 
replaoed  by  their  timeaveraged  values.  The  result  is,  that  this 
approximation  is  justified  from  the  practical  point  of  view 
for  all  machining  operations,  whero  the  borderline  of  stability 
is  low,  that  means  especially  for  those  cases  boing  practioally 
crltioal  and  interesting.  With  regard  to  the  stability  analysis 
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to  bo  performed  quickly  and  economically,  in  the  following  only 
the  time-averaged  values  of  the  direction  coeffi- 

cients  shall  be  considered.  These  values  can  easyly  be  calcu¬ 
lated  by  integration  of  the  direction  functions  between  the 
limits  entrance  angle  and  exit  angle, 

Rxxm  •  at  •  sin  a*  {-  sin  ft  ( cos  2ft  "  cos  2  Tf ) 

+coS/3[2(rA  "fe)-(sin2ft-sin2ft)]} 

Rxyma  -Sjr  ‘  COSot  '  sinx  { -cosft(cos2  fA  -cos  2  ft) 

*  sin  ft  [  2  (ft-  fe;+  (sin 2  ft- sin 2  ft))] 

Rxzm  a  jljf'cosct  ■  cos  x  (-cos  ft  (cos  ft  -  cos  ft ) 
t  sin  ft  (sin  ft  -sinft)  ] 

Ryxm  •  'fif  '  COS  cC  •  Sin  de  {-  COS  ft  (cos  2  fA  -  cos  2  ft) 

*  Sin  ft  C  2  (ft- ft) -( sin  2ft  -  sin  2  (>'9) 

Ryym  *  '  COSdC-  Sin  »  f  sin  ft  ( COS  2ft  -  CoS  2  ft  ) 

+  cos  ft  r.  2  (ft-  ft)  +  (sin  2ft  -  sin  2  ft))} 

Ryxm  «  ~r  cosoC •  cos  af  [sin  ft  ( cos  ft  -  COS  ft  ) 

+  cos  ft  (  Sinft- sinft)  J 

R*xm  K  J~  -sinc(-  Sin  3f  (-  (  Cosft-  coSft)] 

Rzynn  s'jrfjr'  '  SinoC  •  Sin  d€  (ji'nft- sin  ft) 

Rzxm  -  ,4  '  cc,i  *''SinC(-(%-ft) 

On  principle  the  same  relationships  can  be  derived  for  the 
turning  process  as  for  the  milling  process,  if  the  components 
of  the  cutting  force  variation  and  of  the  chip  thickness  variation 
are  defined  as  shown  in  Figure  1,8,  However,  different  simpli¬ 
fications  concerning  the  direction  coefficients  can  be  observed 
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for  the  turning  process,  The  most  important  different)®  la,  that 
they  lire  now  oOnalaiit  valut-5,  bSGRUSS  tt.S  Kr.gls  *f  is  SCKStiir.t . 
Supposing 

f  -  00“ 

which  can  always  bo  satisfied  by  the  reasonable  definition  of 

the  co-ordinate  axes,  It  ,  H  ,  ft  are  aero  Rnd  equations  (1,5) 

a/  i  y  zy 

became  n  .  _  n  ,  .  , 

=  tOSfd  ‘  CCS  «X  '  sm  aP 


Rkl  n  COS/3  ■  C‘dS  o(  1  cor# 

=  -sm  ft  c,osc(  ‘  if/n  06 


Ryx  a  -  r/n/5  ‘  COSO^  •  COS  3C 
c  jinoC'  sm  w* 


R  t  •%  »  ^/no^*cor# 

For  the  orthogonal  turning  process)  givoh  in  figure  1,9,  whore 


only 


I? 


xx 

yx 


0° 

0o° 


-  COB  ft 

“  -sin/9 


f/.U) 


have  to  be  taken  into  aaoount  in  the  chatter  investigation. 


As  already  mentioned  the  transfer  behavior  of  the  cutting 
process  is  described  by  the  so-called  chip  thickness  coeffi¬ 
cient  k0 ,  which  represents  the  rolntionuhip  between  n  chip 
thickness  variation  and  the  corresponding  dynamic  cutting  force, 
Various  efforts  to  determine  this  coefficient  theoretically 
taking  into  account  tool  geometry,  cutting  speed,  feed  rate, 
properties  of  tool  and  workpiece  material  did  not  succeed  until 
now.  However,  the  experimental  duterminaf ion  of  the  k^-values 
is  extremely  difficult  too,  whereby  the  greatest  problem  is 
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to  simulate  na  well  as  possible  the  aoealled  "wave  on  wave  out- 

^  ^  —  —  U  *«<4i  <4  It  A  -  -  •  ««  A  ft*  4  4m*  4  •*  «*  -**  4**  ^  *4  4*  m  m  4  ft-  ft—  «— •  4—  «4-  4ft  4  —*  4  14  -  ■  ,  r,  ■ .  *■»  ^  *j  u  |s|4.  ^  ,  *»  .-, 

VAiijj1  wiuwn  in  «  mh  tiny  x  11  n  uimiiei'iii^  mnuuuuiiD  JJ1  uu  trmn  »  x  nwir- 

fore  the  initial  dynamic  nutting  force  measurement*!  under  thin 
contract  have  been  conducted  for  wave  cutting  and  wave  removing 
only,  which  are  in  general  Bimplifioatibna  of  the  actually 
existing  conditions  (,‘12,133  .  Under  Contract  A?  81(O&2)~900  f  19] 
a  test  rig  to  be  used  on  a  lathe  Was  devetopnd,  by  means  of 
which  also  wave  on  wave  cutting  tests  could  be  performed,  The 
most  important  result  of  the  first  measurements  was,  that 
the  phase  angle  between  the  two  wavy  surface  contours  is  df 
high  influence  upon  the  picked  up  values t  Because  of  this  fact 
the  test  rig  had  to  be  altered  in  such  a  manner,  that  the 
phase  angle  could  ,j e  varied  easyly  and  kept  constant  for  ohuh 
cutting  test  respectively,  Those  alteration!!  required  such 
a  time  delay,  that  no  reliable  and  valid  results  could  be  ob¬ 
tained  until  now, 


The  applicability  of  the  random  hoise  theory  in  dynamic  dotting 
teats  proved  alec  impossible  until  now  due  to  the  effect  of 
the  phnae  angle  fl9}  , 

However,  from  the  outtihg  tests  carried  out  up  to  this  time 
C 12 i id , 10 i 0l  some  important  results  could  be  obtained,  The  most 
interesting  fact  is,  that  the  generally  Used  equation 

k0  "  f  (1.10) 

where  the  outtihg  force  variation  is  regarded  ae  directly  pro¬ 
portional  to  the  chip  thickness  variation,  is  an  approximation, 
which  Is  permissible  under  aurtain  conditions  only,  This 
equation  is  mainly  valid  for  relatively  low  chatter  frequencies, 
For  higher  frequonuimi  it  must  bo  considered,  that  k0  is  depen¬ 
dent  on  the  frequency  with  respect  to  amplitude  and  phase  angle, 
Thus  follows 
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Equation  (1.11)  related  to  the  width  of  out  w  yields 

kcw(Lui)m  A  (iw)  B  JL.  t(.u)  (/./l) 

Herein  the  related  chip  thickness  coefficient  k^Ci  w  )  reprs** 

0  •* 

eente  a  special  property  of  material ,  characteristic  for  the 
chosen  cutting  conditions.  As  already  pointed  out,  the  Invest!** 
gat  ions  performed  in  this  program,  are  pritaaryly  directed  at  thu 
heavy  duty  machining  of  the  usual  Fe*-matertals  with  nutting 
speeds  in  the  range  of  about  70  , , ,  200  m/ain.  Due  to  this  the 
so-called  low  speed  stability,  which  is  preferably  interesting 
in  the  machining  of  high  strength  and  hard  materials  C 9]  ,  can 
be  let  out  of  consideration.  Over  the  cited  speed  range  the 
chip  thickness  coefficient  is  therefore  taken  as  independent 
from  the  cutting  speed. 


The  chip  thickness  coefficients  determined  in  an  orthogonal 
plunge  turning  teat,  can  bs  applied  to  face  Hilling  according 
to  ft  If  D»*a  1*9.  Fro«  several  results  of  dynamic  cutting  tests 
it  can  be  seen,  that  the  spatial  components  of  the  cutting  force 
*  especially  at  higher  frequencies  «  have  distinguished  phase 
shifting  to  the  chip  thickness  variation,  this  means,  that  the 
direction  of  the  resultant  dynamic  cutting  force  is  time -Variant 
in  space i  With 
I'll  l*  main  cutting  force  and 

-  thrust  force  consisting  of  a  radial  and  an  axial  component, 
the  transfer  behavior  of  the  cutting  process  can  bs  described 
by  the  equations 


*kH 
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and  related  to  the  width  of  out 
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Aooordihg  to  the  equations  <  1 1 1 3 >  and  (1,14)  itt  Figure  1,10  the 

transfer  behavior  of  the  cutting  process  is  represented  by  each 
two  blocks  -  phasefree  blocks  w»  kQWH(iu  )»  W* 
and  phase  shifting  blooke 

(u)  B  (cosVy  (uj)+  i  jih  Vw  ltd)) 

tuj)  «  (  cos  lui)  -t  i  jbr>  lid)) 

for  the  two  oomponents  Fjj  and  FA.  Projecting  these  components 
into  the  direction  of  the  resultant  cutting  force  by  the  an¬ 
gular  functions 

R^(m)  "  A  m  i9  ( u*)  cosc((u' )  (>•  IS  ) 

Ra  aol  «  sfhK  ‘caj3(vj)  1  t'Oiti  (oj!  ('  IJi*  ) 

which  can  be  derived  from  Figure  1,0,  and  summation  yield  the 
cutting  force  F(tuj),  whiah  reacts  upon  the  vibration  system 
of  the  machine,  As  it  man  bs  seen  on  Figure  1,10,  in  an  aotual 
out  ting  process  the  width  of  out  w  ia  the  gain  factor,  by 
which  the  amount  of  the  dynamic  cutting  force  and  therewith  the 
stability  id  direotly  influenced .  Tho  definitions  of  w  in  tur¬ 
ning  and  milling  is  given  by,  Figure  1,11, 

If  tho  direction  of  the  dynamic  cutting  force  is  timevariant, 
the  ratio  of  the  meaaured  cutting  force  components  and  there¬ 
with  the  angle  $  and  all  values  connected  with/9  are  also  time- 
variant,  Taking  into  account  this  time  variant  behavior  in  the 
stability  analysis  would  make  the  calculation  difficult  and 
very  time  consumptive.  With  rogard  to  an  economical  performance 
of  the  chatter  investigation  it  shall  therefore  bo  assumed  in 
the  following,  that  tho  cutting  force  direction  may  be  depen¬ 
dent  from  the  frequency  but  time  invariant  for  each  discrete 
frequency.  This  assumption  may  uause  slightly  incorrect  results) 
however,  as  demonstrated  ir,  the  chapters  1,3,1  and  the  ob. 
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triiiSu  results  are  satisf&sterily  a  as  at  far  the  practical  use. 
The  chip  thickness  onef riclents  used  in  the  following  calcula¬ 
tions  ,  show  by  far  time-invariant  aufcting  force  direct  ions 
Within  the  interesting  frequency  range  [la  ,13,0]  ,  Here  the 
maximum  phase  angles  between  the  force  components  come  up  to 
about  10°,  This  result  eorrolntcs  quite  well  to  the  results  of 
dynamic  cutting  tests  conducted  by  Poincek,  of  which  some  data 
were  discussed  at  the  0,1  ,n,P, ‘-Meeting  of  Working  Group  "Ma" 
oh  1  February  i960  in  Paris  (France),  Figure  1,12  shows  a 
schematics!  drawing  of  an  orthogonal.  pluhge  turning  process. 

Due  to  the  chip  thickness  u(iin)  the  static  cutting  force  F0 
Is  superimposed  by  a  dynnwia  cutting  force  F(t),  of  whioh  the 
direatioh  must  not  correlate  to  the  direction  of  FQ,  In  the 
given  example  the  ratio  of  the  foroe  oomponents  shall  be 
1‘1(/F^  -  a.  In  the  cane  "  V'H  "  °°’  monnB  no  Phase  shif¬ 

ting  between  the  components,  tho  direction  of  the  cutting  force 
and  thus  the  angle  0  are  constant,  it  follows 


fi'circto^L  (1.17) 

Is  there  uhy  phase  difference  of  for  instance  ■  10° 

(interrupted  curve  of  Fj{),  the  end  of  the  dynamic  cutting  force 
Vector  follows  an  elliptical  curve.  However,  as  it  can  be  seen 
from  the  drawing,  for  a  small  phase  difference  tho  deviation 
of  the  cutting  foroe  direction  from  the  time  invariant  direc¬ 
tion  is  only  small,  When  the  outting  foroe  veotor  itsulf  is 
great ,  This  moans,  that  in  the  oases  of  small  phase  differences 
oqu,  (1,17)  is  valid  with  good  approximation.  For  toolH  with 
■at  -  DO*’  one  ohtaines  according  to  Figure  1.9 
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Cut*  and  Several  Degrees  of  Freedom 


In  connection  with  tho  derivation  of  the  dii'eotion  coefficients 
the  relationship  between  a  relative  displacement  at  the  cutting 
point  and  the  corresponding  cutting  force  was  detailed,  In  ma¬ 
chining  processes  this  force  reacts  upon  the  machine  structure  * 
so  that  the  process  can  he  represented  by  the  closed  control 
loop  in  Ftguro  1,13,  if  a  vibration  system  with  one  degree  of 
freedom  and  real  behavior  of  the  cutting  process  arc  supposed, 
Hereby  the  transfer  behavior  of  the  cutting  process  is  represen- 
tod  by  the  material  characteristic  value  kRW  and  the  gain  factor, 
the  deapth  of  out  \v,  Assuming  real  behavior  of  the  oUttihg  pro¬ 
cess,  equ.  (1,10)  yields 


For  a  system  with  three  degrees  of  freedom  the  block-diagram 
in  Figure  1,14  is  obtained.  Here  the  total  chip  thickness 
variation  u^  at  the  i-th  tooth  results  from  the  components 
U^(x),  tt^(y),  (z) i  which  are  caused  by  the  displacements 

x,  y,  7,,  Hy  (w  •  hQW)  the  chip  thickness  variation  produces  a 
cutting  force  variation  F^ ,  The  components  F^,  F^yl  F^  of  F^ 
act  upon  the  machine  structure  again,  ny  reasonably  combining 
the  blocks  of  tho  control  loop  and  considering  (1,0)  the  dia¬ 
gram  in  Figure  1 .15  is  obtained, 

With 


Ac 


H, 


+  Py„  A 


and 

Ko 


+  R  72  •  GtX 


(1.21> 

(1.22) 


tho  very  simple  control  loop  in  Figure  1,10  results, 

Horo  tho  index  0  indicates,  that  tho  additional  phase  shifting 
duo  to  tho  complex  behavior  of  the  cutting  process  is  not  yet 
taken  into  account.  The  influence  of  this  complex  behavior  is 
dotal  led  in  chapter  1,2,3, 
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1.2  Sfcabil i fry  Analysts  by  Means  of  tho  Nyqulgtt.-'lrltcvlon 


As  it  was  found  out,  tho  dynamic  maohinlnrc  process  onn  always 
-  oven  for  systems  with  several  dorrroos  of  freedom  -  he  reduced 
to  a  nonintermoshod  closed  control  loop,  This  fact  it  very  inter¬ 
esting  and  of  high  Importance  for  t bo  stability  analysis,  As  tho 
transfer  behavior  of  tho  machine  and  the  cutting  process  is 
available  in  the  form  of  experimental ly  measured  response  loci 
only,  the  stability  analysis  is  suitably  performed  applying  tho 
well-known  NyqUist-Criterion ,  fty  means  of  this  criterion  the 
stability  of  the  closed  loop  is  analysed  bn  the  basis  of  tho 
response  locus  F  (iw)  of  tho  open  loop,  (riven  in  Figure  1,17. 

1,2,1  Control  L bop  Eq untion 


With  z  input. 

O 

Z  "  output 

the  equation  qf  tho  open  loop  response  locus  is  obtained 


to  (iu) 


(1 .23) 


For  tho  Chatter  atinlysis  the  total  response  loous  must  not  be 
known ,  because  it  is  only  important  to  determine  those  condi¬ 
tions,  whore  the  boi’der  of  stability  is  readied,  At  the  stabi¬ 
lity  border  (She  output  «l  is  identical  l:o  the  input  /,Q  concerning 
phase  And  amplitude: 

»  1  (1.24) 

/jo 

Thus  the  orosspoint  of  the  response  locus  F  with  the  real  axis 
is  tho  point  (+  1,0),  the  so-called  critical  point.  As  to  the 
chatter  this  moans,  that  vibrations  once  excited  donl;  neither 
inoroftwo  nor  decrease.  With  a  orosapoint  below  (+  1)  the  system 
is  nbsolutjy  stable ,  a  orosspoint  above  (1,0)  means  instability 
of  the  system,  Equ,  (1.22)  and  (1.21)  ylol.d  the  omul  It.  ion  for  the 
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border  of  stability: 

*o  “  Ro  u  ' 

1.2,2  Detsrmination  of  the  Stability  Borderline 

By  mean*  of  this  equation  for  eaoh  frequency  the  value  Kc  oan 
bw  calculated,  for  which  the  control  loop  reaches  the  stability 
border,  The  value  K0  Includes  the  ohlp  thickness  coefficient, 
on  which  the  deapth  of  out  Is  immediately  dependent  (1,22  and 
1.20),  as  well  as  the  delay  time  T  ,  which  is  used  for  deter¬ 
mining  the  corresponding  rotational  speeds  of  the  workpiece 
(in  turning)  and  of  the  cutter  (in  milling): 


where  T  "  time  for  one  revolution  of  the  workpiece 

(turning)  and  the  cutter  respectively  (milling) 

Z  -  total  number  of  teeth  on  the  autter 


The  values  AQ  and  KQ  are  oomplex  magnitudes,  which  can  be 
written  as 
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(A0r 

+  i 

Ai  ) 

(1.27) 

K(>  * 

(  Kor 

*  i 

Koi ) 

(1.28) 

with  Aof 

“  M0} 

A{,  i  ■  Jm  {  A0) 

Kor 

*  ft  (  Ko} 

Koi  ■  jm  f  Kg] 

Equ.  (1.28) 

compared  1 

:o  (1 

,22 

)  yields 
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+  i  Koi  )  ■ 
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kc»  ‘  (  A) 

(1.29) 

and  with 
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cos  cj T  -  i  j. in  wt 
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(Kar  +  I  K«i ) 


( e©<  u  T  -  A  ~  i  sin  w  T ) 


(hV>) 


By  this  equation  (w  ■  kflw)  and  T  can  be  determined  for  a  Riven 
frequency  w  .  For  this  purpose  equ.  (1,30)  is  expanded  and 
arranged  to 


Kpr  "  i  Kqi 
(K‘  +  Kii) 


tv  ■  kCw  ( cos  cj T  “  A  -  i  sir i  yT) 


Prom  this  equation  the  real-  and  imaginary  parts  of  Kc  are 
obtained : 

Kor  *>  JV.  kcw  (  Kor  +  Koi )  ( con  u>T -  <1 ) 

Koi  -  W-kcw  (K*f  +  «oi  )  (  sin  tot) 

(tV'kcw)  assumed  to  be  real,  as  mentioned  above,  yields 
by  division 

J?±L  -  'sinuT  Uu) 

Kor  CotoT-'l 
and  with 

c to  ULLs 

3  i  T~  cos  tot 


and  considering  (1.23) 

w  T  »  2  arc  ctg  (-  -j5— ) 


(I.  32) 


This  equation  describes  an  inverse  trigonometric  function,  of 
which  the  main  values  are  limited  by 

0  -  to  T0  -  P  7T 

Considering  the  ambiguity  of  equ.  (1.32)  the  general  solution 
can  be  formulated  as 


U)  T  -  <*>  T_  +  ?J  m 


(nr)  z  0,  1,  2,3  “•  ) 


r*T*  ''T-  ^ 

T'  T»  ♦  — ~ 


(1.33) 
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with 


M. 

7  e 

the  well  known  relation 


(1.34) 


is  obtained,  by  means  of  which  the  rotational  speed  correspon¬ 
ding  to  the  frequencies  in  the  response  locus  oan  be  calculated, 


Corresponding  to  the  pair  of  values  T  and  (w  •  ko^)  whiah  satis¬ 
fies  equ,  (1.22),  the  value  (w  1  kQW)  can  be  determined. 

Equ.  (1.22)  expanded  by  K0|>  yields  after  some  arranging 


tv  •  ktw  •  Ktf 


( 4  + 


( cos  o>T‘ 

Kor 


4  -  i  Ain  a)T ) 


and  with  (1.31)  after  suitable  expanding  and  arranging  the 
relation 


W  •  AC(y 


(1.38) 


is  obtained.  As  detailed  above,  the  system  oan  only  beoame 
unstable,  if  the  response  loous  of  the  open  loop  crosses  the 
real  axis  at  a  point  >  (+  i  ),  Considering  (1.35)  it  follows 
from  this  fact,  that  Kor  must  be  negative,  because  only  a  posi¬ 
tive  ohlp  thickness  coefficient  and  a  positive  deapth  of  out 
respectively  have  a  practically  reasonable  meaning.  Otherwise 
the  system  is  stable. 


With  (1.2$  and  (1.28)  equ.  (1.25)  becomes 

(  Aor  +  t  Aoi)  -  (  Kor  ■+  i  ■  Hoi  )  -  0  (1.36) 

This  means 

V  -  Kor  "  0  <1-30a> 

A0l  -  Kol  -  0  (1.36b) 

Considering  (1.21)  it  follows  from  (1.36a) 

K0Ka  Re{RMK'Gxx+*,'  +  •  Gj }  (1.37) 


-  23  - 


and  with  (1.35) 


IV  a 


,1 _ 

2  kew  •  Hg  [  R„x  •  6km  +  •  •  •  +  R/7*6?2] 


(1.38) 


the  relation  for  the  calculation  of  the  aotual  critical  deapth 
of  out  results.  As  it  can  be  seen,  therewith  the  stability  ana¬ 
lysis  of  machining  operations  can  in  every  case  be  reduced 
to  an  equation,  which  is  easily  to  be  solved  by  hand. 

.1,2,3  Complex  Transfer  Function  of  the  Cutting  Process 


As  detailed  in  chapter  1.1.2,  with  regard  to  easy  performance  of 
the  stability  analysis  the  direction  of  the  cutting  force  in 
space  is  assumed  to  bo  time-invariant,  that  means 

(oj)  ■  ^  (tu)  =>  H’ i'ui) 

Thus  the  complex  transfer  behavior  of  the  cutting  process  oan 
be  represented  by  the  block  diagram  In  Figure  1,18  a.  From  the 
diagram  the  relationship  between  the  chip  thickness  variation 
u  and  the  cutting  force  variation  F  oan  directly  be  derived: 


Lduj)  z  b/\l<ch,H\(oj)-  + 


With 


R  =  C  Rh(Oj)  +  (w)1  ra  (u,j1 


^A\  <u)  ■  Ra  (to)]  1  ( CoS  V (U>)4-  i  J fin  ^(u)'j 

(1.19) 


the  diagram  is  simplified  to  the  one  on  Figure  1.18  b.  As  it 
oan  be  seen  here,  the  transfer  funotion  of  the  cutting  process 
includes  proportional  factors  w,UowjiI  (w)»  anda  phase 

shifting  block  (cos  V(u|+  i  sinYluJ)  ).  Considering  this  be¬ 
havior  in  the  block  diagram  of  the  machining  operation  yields 
Figure  1,19  a.  If  the  phase  shifting  blook  is  combined  to  the 
block  AQ  -  as  shown  in  Figure  1,19  b  -  it  follows 
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A  Uoj)  •  A0  (iwj  •  (toi  “'fCot)  +  Cjnh  V  (u)) 


(1.40) 


This  equation  oan  be  Interpreted  In  such  a  way,  that  each  vec¬ 
tor  of  the  machine  response  loci  rotates  against  the  coordinate 
system  by  that  phase  angle  ,  which  is  valid  at  the  concerned 
frequency  wi  .  In  the  machine  response  loci  0(1  cj)  the  angle  (■' ft ) 
means  a  phase  lag  of  the  displacement  against  the  force.  Now 
an  angle  (-Yi;  )  increases  this  phase  lag,  as  shown  in  Fig,  1.20. 
The  other  blooks  of  the  loop  oan  be  written  in  shorthand  nota¬ 
tion,  analogly  to  (1.22) 

R  •  kcvbt  *  (e~lwr-yi)  (1.41) 

so  that  the  control  loop  in  Figure  1,19  c  results. 

The  stability  analysis  is  now  carried  out  in  the  same  manner 
as  shown  above  for  the  case  of  a  real  chip  thickness  coeffi¬ 
cient,  With  the  control  loop  equation 

A  -  K  -  0  (1.42) 

the  equation  for  calculating  the  oritioal  deapth  of  out  is 
obtained: 


2  •  R  •  kcwH  '  (  fix*  ’  +  "  •  ■+*  }  (1.43) 
The  rotational  speeds  of  the  workpiece  and  the  cutter  respec¬ 
tively  can  be  determined  by  means  of  equ,  (1.34)  again. 

1.3  Practical  Chatter  Investigations 


By  the  above  described  relationships  it  is  now  possible  to 
perform  the  theoretical  stability  analysis,  if  the  charac¬ 
teristic  values  and  the  transfer  behaviour  of  the  cutting 
process  and  the  machine  are  known.  As  already  mentioned, 
the  cutting  process  data  are  only  dependent  on  the  cutting 
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conditions;  thus  they  cnn  bo  determined  In  separate  cutting 
fore®  measurements  t\nd  are  than  available  for  further  sta¬ 
bility  investigations,  Against  that  the  machine  response  loci 
arc  characteristic  data,  which  arc  in  (tenoral  specific  for 
each  machine.  Because  of  this  the  response  loci  must  be  picked 
up  at  each  machine  to  bo  investigated  with  regard  to  machining 
stability . 

During  the  time  covered  by  this  contract,  17  lathes  of  diffe¬ 
rent  typos  and  manufacturers  and  altogether  0  milling  machines 
-  of  which  two  medium  sized  vertical  milling  machines,  four 
single  column  boring  and  milling  machines  and  three  large  portal 
milling  maohines  -  have  extensively  been  investigated,  The  most 
important  results  of  those  investigations  have  been  published 
in  the  Quarterly  Technical  Reports  under  this  contract  (see 
[1.1  -  18!)  ),  The  general  purpose  of  these  machine  investiga¬ 
tions  was  to  verify  the  validity  of  the  theoretical  research 
work  by  means  of  experimental  machining  tests  on  the  one  hand 
and  on  the  other  hand  to  find  out,  how  the  chatter  analysis 
can  be  performed  under  consideration  of  the  different  machin¬ 
ing  problems  existing  in  prnotico.  Some  representative  investi¬ 
gations  and  results  shall  be  detailed  in  the  following, 

At  first  it  will  be  shown,  how  the  stability  borderline  aan 
easely  be  determined  by  hand  especially  for  turning  processes, 
Hereby  the  performance  of  the  calculation  ann  considerably 
bo  simplified  by  means  of  graphical  illustrations.  For  exten¬ 
sive  and  complete  chatter  investigations  in  general  the  follow¬ 
ing  parameters  have  to  bo  taken  into  account: 

a)  concerning  the  machine: 

working  range 

range  of  rotational  speed 

installed  power 

maximum  torsional  moment 
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b)  concerning  the  machining  process; 

workpiece  material  -  tool  material 

milling  breadth 

diameter  of  tho  outter 

number  of  tooth  on  the  tool 

tool  geometry 

configuration  workpiece  -  tool 

Considering  theae  numerous  parameters  the  chatter*  investigation 
cannot'  be  performed  with  economic  calculation  effort  without 
using  electronical  computers,  Therefore  with  regard  to  tho 
rational  calculation  of  stability  charts  digital  computer 
programmed  have  boon  developed,  tho  application  of  Which  in  do- 
monstrntod  by  an  cjcamplo  of  ohnttor  investigations  on  rv  vurtlcal 
milling  machino. 

1.3.1  Chatter  Investigations  on  a  Latlio 

These  chatter  investigations,  Which  included  tho  dynamio  monsuro 
ments,  tho  calculation  of  the  critical  width  of  out  us  well  as 
for  comparison  experimental  cutting  tests,  have  been  conducted 
at  a  modium-aissed  latho  (30  Ml),  The  calculations  shall  be  li¬ 
mited  to  the  orthogonal  plunge  cut,  because  hereby  very  easily 
understandable  geometrical  conditions  are  existing. 

1 . 3 , 1 ,  .1.  Dynamio  Behaviour  of  the  Investigated  Lathe 

For  the  dynamic  investigation  of  the  machine  the  contactless 
electro-magnetic  oxciter  was  applied,  which  has  been  developed 
under  this  contract  (see  [17]  ).  The  required  response  loci, 
picked  up  for  a  rotational  speed  of  the  spindle  of  n  *>  140  min’"1 
are  plotted  in  Figure  1,31.  It  can  bo  stated,  that  the  machine 
has  two  main  vibration  modes!  a  bending  mode  of  the  spindle 
at  lie*  und  a  rocking  mode  of  tho  tool  head  at  170  cps. 
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'n?  tilts  wtnbiUty  nnalyuin  tt  In  r*  f  litirh  4  m  tint1  ban  no . 
that  the  compliance  in  the  bunding  mode  at  .1.16  ops  with  the  ro¬ 
tating  spindle  overcomes  the  value  for  the  stand  a till  spindlG 
by  about  1.00  %.  Tho  modo  at  170  ops  is  naturally  not  in- 
fluonood  by  the  spindle  rotation, 

1.3,  1,2  Determination  of  the  Gritioal  Width  of  Cut 


When  manually  performing  tho  stability  analysis,  tho  calcu¬ 
lation  uffort  is  considerably  dependent  on  tho  transfer  be¬ 
haviour  of  tho  cutting  praaeas.  If  there  is  a  frequency-de¬ 
pendant  phase  shift  existing  betweon  tho  chip  thiaknoss  - 
and  the  cutting  foroo  variation,  then  each  vector  in  the 
rosponso  loci  has  to  be  turned  by  tho  corresponding  angio  V** 

As  show  in  Figure  1.20  a  negative  angle  Y  means  a  phasolag 
of  tho  cutting  forao  -  against  the  chip  thickness  variation, 
i.o.  tho  eomplioanco  vectors  have  to  bo  turned  clockwise 
against  the  coordinate  system, 

This  procedure  oil  principle  being  simple  but  in  prnotioe  being 
relatively  time  aonatuuptiTO,  aan  considerably  be  simplified 
if  the  phase  shift  is  oonstant  within  the  interesting  fre¬ 
quency  range,  In  suoh  oases  the  phase  angle  is  idontiaal 
far  all  voators  of  the  rosponso  loci,  so  that  only  the  coor¬ 
dinate  system  must  be  turned  -  counter  clockwise  forV‘<  0°.  The 
chip  thiakness  coefficient  used  in  the  following  example,  has 
only  negligible  phase  shift  in  the  interesting  frequency  ranee. 

It  was  determined  in  [1.3]  , 

kcw  "  3»°  *  10^  kp/cm3 

ft  «  67^ 

Thus  equ.  (1.38)  can  bo  used  for  tho  determination  of  the  critical 
width  of  cut.  As  already  mentioned,  chatter  can  only  occur,  If 

Re  /Rxx  ’  °xx  +  Ryx  ‘  Gxy}  <  0  ftt  one  or  more  fl,eclue,,cieH* 
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if  this  is  sot  satisfied,  th®  gy»t*m  wi.I1  h»  stable  for  what¬ 
ever  high  widths  of  cut.  The  above  summation  is  reasonably  be  per¬ 
formed  graphically.  Figure  1,22  shows  the  run  of  the  functions 
Ho  {  Rxx  and  Re  £Ryx  _  ^xy  3  as  woi:l  ftl9  thG  resultant 

curve  versus  tho  frequency.  The  resultant  curve  shows  three 
minima  with  negative  sign,  tho  absolutoly  greatest  of  which 
with  (  “0,2)  yum/kp  occurs  at  123  ops,  For  this  the  minimum 
critical  width  of  out  is  calculated} 

w  «  0,4  nn 

For  tho  two  other  minima  at  160  cps  and  200  cps  widths  of  out  of 
21mm  and  24  mm  respectively  are  obtained,  which  however  are  not 
interesting  from  the  praotioal  point  of  view.  The  borderlines 
of  stability  are  plotted  in  Figure  1,23.  As  it  can  bo  seen, 
tho  theoretical  value  of  0,4  mm  is  quite  well  verified  by  the 
additionally  plotted  experimental  results: 

If  tho  stability  calculation  would  be  performed  using  the  equa¬ 
tions  (1.36)  and  (1,34),  the  theoretical  stability  chart  would 
show  the  aharaateristloal  unstable  lobes,  the  so  called  "chatter 
lobes."  The  minimum  of  these  lobes  is  Identical  to  the  above 
calculated  critical  width  of  cut  of  8,4  nun.  However  in  turning 
operations  these  lobes  overlap  so  far  in  general,  that  the  ro¬ 
tational  speeds,  where  \  the  minimum  width  of  cut  is  reached 
lie  very  closely  to  one  another  in  the  usual  range  of  rota¬ 
tional  speed.  Thus  in  praotioe  the  width  of  cut  to  be  reached 
is  equal  to  the  minimum  critical  width  of  out.  For  this  reason 
the  calculation  of  the  -  theoretically  -  critical  rotational 
speeds  can  be  given  up  in  general  for  turning  operations. 

1.3.2  Chatter  Investigations  on  a  Milling  Machine  Using 
tho  Digital  Computer 

1 . 3. 2. 1  Dynamic  Behaviour  of  the  Investigated  Milling  Machine 


The  investigations  were  carried  out  on  a  modium-Bized  vertical 
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milling  machine  (  llltVf).  Tho  determination  of  the  response  leal 
required  for  tho  stability  analysis,  was  performed  by  convon- 
tional  exciter  testa  nt  tho  standing  still  machine.  In  Figure  l.iM 
the  response  IogI  for  tho  following  working  position  aro  plotted: 

quill:  out 

spindle :  above  the  tablo  centre 
cutter:  70  mm  above  tho  tablo 

As  it  can  be  soon  in  the  figure  tho  machine  shows  mainly  four 
vibration  modes  in  the  tested  frequency  range,  namely  at 
B7  ops,  G3  ops,  68  ops  and  138  ops.  These  aro  illustrated  in 
Figaro  1.25.  The  vibration  modes  oxoited  in  y-direation  at 
65  ops  and  138  cps  aro  also  excited  by  F„  at  tho  same  fro- 

Cl 

quencies.  As  it  was  stated,  tho  vibration  modes  depend  con¬ 
siderably  on  the  position  of  the  travelling  machine  elements 
concerning  the  cowplianoe  and  tho  natural  frequency.  Dy  the  po¬ 
sition  of  the  table  and  tho  cross  slide  respectively  in  the 
xy-plane  oxpocially  the  vibration  modes  at  57  ops  and  65  ops  are 
influenced.  As  shown  in  the  following,  preferably  these  modes 
are  decisive  for  tho  chatter  behaviour.  The  variation  of  tho 
machine  dynamics  duo  to  the  position  of  the  elements  has  ex¬ 
tensively  been  discuosod  in  [17], 

1 . 3 . 2 , 2  Development  and  Application  of  the  Computer  Programmes 
HATTER  1.  and  RATTER  2 


Corresponding  to  the  thaoretical  derivations  in  chapter  1.2  the 
calculation  of  a  stability  ohart  includes  the  following  stops: 

1.  Calculation  of  the  outting  force  angles  d  (icu)  and/?(iw); 

2.  Calculation  of  time-invariant  direction  coefficients 
U  ...  \\zz  and  of  the  coefficients  R(r ,  R^; 

3.  Calculation  of  tho  products  II  •  G  (iu)  ...  g  .  Gr,„  (iu) 

aa  2315 

and  summation  at  each  frcquoncy;  the  results  can  bo  plotted 
ns  tho  response  locus  Aq  (lu); 
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4.  Multiplication  of  AQ  (iu>)  by  the  phase  shifting  block 
"H1  (leu);  this  yields  the  response  locus  A  (iCJ); 

5.  Calculation  of  the  chip  thickness  coefficient  k0jj  and 
of  tho  critical  deapth  of  cut  at  all  f requencics , 
where  A  (iu>)  has  a  negative  real  part} 

G,  Calculation  of  the  rotational  apaoda  of  tho  cutter 
within  an  interesting  range  of  the  cutting  speed} 

7 4  Plot  the  stability  chart  and  /  or  store  the  results, 

Hy  moans  of  tho  Computer  program  RATl'DR.  1,  tho  flow  chart 
of  which  is  given  in  Figure  1 ,2(3,  the  above  listed  stops 
can  be  carried  out.  In  tho  given  form  by  thin  program  one 
complete  stability  chart  fox*  tho  pai'ninetox'd  in  tile  hoftd  of 
the  program  is  obtained  by  each  calculation  run.  The  in- 
veutignlionu  Uotnilod  hero,  wore  directed  to  the 

left  end  position 
middle  position 
right  and  position 

of  tlui  table  for  constant  positions  o;F  the  quill  and  tho  knee. 

In  these  three  positions  stability  charts  have  boon  calcula¬ 
ted  for  the  four  milling  ootifigurutlonts  shown  in  Flguro  ,i ,  37. 

As  it  can  bo  seen,  those  four  configu, rut  ions  roprosont  four 
different  possibilities  to  perform  ono  given  milling  oporntion 
at  the  workpiece. 

In  Figure  tjJjJJ  tho  stability  uhni’ts  for  tho  four  oonf igurntions 
in  the  throo  table  positions  nvo  plotted,  Comparing  tho  sta¬ 
bility  charts  it  can  be  stated,  that  considerable  ditfovonuon 
in  the  critical  deapth  of  uxit  occur  on  tho  one  hand  for  different 
configurations  in  one  tablo  position  and  on  the  other  hand 
also  for  tho  snme  configuration  in  different  tabid  positions, 

The  last  comparison  is  practically  in to run ting,  when  tho  given 
machining  operation  requires  n  long  travel,  of  the  tablo,  Hero 
the  minimum  stability  is  thou  decisive  for  tho  maximally 


31 


po.MfjUilo  doapth  of  cut .  Foe  inulrvnco  with  an  entranoo  angle 
of  i uu°  n  oiinttoriree  deapth  or  out  of  only  >,0  mm  can  bo 
reached,  although  0,[i  mm  ahattarfroo  doapth  of  out  Is  poauib- 
lc  In  the  loft  end  pout  Lion.  Foe  all  table  positions  tho  moot 
favourable  conf iguration  id  given  foe  an  ontranao  angle  of 
90°,  because  chatter  must  bo  cxapocted  lioro  not  bo  low  t),0  mm, 
lit  the  theoretical  calculations  n  chip  thiohness  ooG'f fioiont  de¬ 
termined  in  j 13]  Who  Utfed t 

RcW  "*  S,0  *  to4  ltp/em2 

ft  0?o 

Oi-  7  0^0 

This  coefficient  shows  only  real  behaviour  ovor  tho  interesting 
frequency  range,, 

For  tho  verification  of  tho  thoorci  licml  results  in  Home  oases 
experimental  chatter  tests  have  boon  carried  out,  the  results 
of  which  are  additionally  plotted  in  tho  ratability  chart,  it 
utth  b<a  istated*  that  thn  correlation  between  tho  theoretical 
and  ttJcjporlmental  results  id  ah  to  the  deapth  of  out  aad 
ha  to  thn  rotational  speeds  and  this  chatter  frequencies  re¬ 
spectively. 

In  order  to  make  obvious,  in  how  far  this  installed  power  of 
the  machine  la  actually  turned  to  profit  and  therewith 
to  get  it  ori  tor  ion  for  tho  ooonomy  of  machining,  it  is 
reasonable  to  plat  the  limitation  of  the  deapth  of  out 
due  to  tho  installed  power  and  duo  to  the  maximum  torsional 
moment  additionally  in  tho  (liability  uhavts.  Tho  doapth  of 
out  limited  by  the  installed  power  can  bo  calculated  for  a 
given  feed  rate  per  tooth  by  tihu  liquation 

_ 
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where  NQ  «  installed  power 

^  -  total  efficiency 


k  .  .  ,n  specific  cutting  forco 

s  V4 

h  ,u  averaged  chip  thickness 

h  «  2.£2_  ~3  1  Jit  jinn  .. 

m  TT  a  (Yt-rg) 

1-53  u  cutting  forcooxponont 

a0  »  number  of  teeth  being  in  contact 

»  -  2  •  <A  -TV° 


From  tabulns 

(  1.  -  sd  )  -  0,88 

2 

k  ^  1BG  kp/mm 

ato  taken. 


Tha  offioiendy  ij  was  determined  In  the  tests) 
^  ■*  0,55 


In  the  range  of  low  spindle  speeds  the  limitation  is  not  given 
by  the  power  but  by  the  maximally  admissible  torsional  moment. 
Assuming  that  the  tangential  force  on  the  cutter  is  quite  in¬ 
dependent  on  the  cutting  speed*  the  doapth  of  out  limitation 
duo  to  the  torsional  momont  ann  bo  calculated  by 

(1,45) 

The  calculation  of  thoso  limitations  can  also  be  performed 
by  the  oomputor  program  KATTUK  1.  an  shown  in  tho  flow  chart, 

IJy  moans  of  tho  formulas 
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mined  lor  whatever  feed  rates. 

As  already  mentioned  tho  program  RATTER  1.  is  able  to  produce 
one  stability  chart  for  a  given  parameter  combination  by  each 
run.  If  more  than  one  combination  are  interesting,  the  program 
can  easily  be  completed,  so  that  tho  corresponding  parameters 
can  be  varied  automatically,  Tho  results  of  such  an  extensive 
stability  analysis  in  certain  circumstances  is  a  groat  number 
of  stability  charts,  which  characterize  tho  chatter ] behaviour 
of  the  machine  under  the  interesting  conditions.  For  the  or¬ 
ganization  of  work  such  datn  may  bo  very  useful,  howover  as  to 
the  application  in  the  workshop  this  groat  number  of  stability 
charts  is  only  difficultly  to  survey  and  to  handle, 

Thus  it  raises  the  problem  to  reduce  the  number  of  data,  but 
without  leaving  out  and  neglecting  the  praatioally  required 
information,  ThlB  aould  be  realized  by  altering  and  further 
developing  the  program  RATTER  1,  as  detailed  in  the  following. 

The  program  RATTER  2  (Fig,  1,29)  was  established  on  the  basis  of  tho 
problem  being  often  existing  in  practice,  that  for  a  given 
machining  operation  at  the  Workpiece  the  optimum  conditions 
must  bo  determined.  Apart  from  special  oases,  as  for  in¬ 
stance  maohining  with  boring  bars,  rough  machining  operations 
are  in  general  critical  with  regard  to  the  chatter  behaviour. 

Here  one  reduction  is  already  possible  dus  to  the  fact,  that 
for  given  workpiece  materials  certain  toolH  havo  been  proved, 

Thus  for  turning  tools  and  outters  tho  tool  matorinl  and  tho 
tool  geomotry  are  given,  Dy  the  workpiooe  -  tool  combination 
hIho  tho  cutting  apoed  range  and  thu  food  rate  per  tooth  are 
limited,  If  a  certain  maohina  and  a  tool  are  chosen  lor  thu 
given  machining  operation,  at  last  thu  most  favourable  po¬ 
sition  in  the  working  range  iim  woll  as  tho  optimum  configu¬ 
ration  workjiloce  -  tool  must  bo  detorminod  by  moann  of  tho 
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Furthor  considerable  reduction  of  Lhu  data  will  bG  possible 

duo  to  the  fact,  that  it  is  preferably  of  interest  in  practice,  i 

which  lowest  chnttorfrec  denpth  of  cut  can  surely  be  realized,  \ 

This  question  can  be  answered  by  specifying  only  ono  value 
nnmoly  the  minimum  doapth  of  cut  in  tho  stability  chart. 

In  Flguro  1,30  some  rasults  of  the  stability  calculation  ! 

by  moans  of  the  computer  program  HATTER  2  are  given.  Here  the  j 

run  of  tho  minimum  stability  of  all  possible  milling  oon-  1 

figurations  for  a  given  n/D-rntio  and  a  cuttor  with  given 

tool  goomotry  is  plotted  against  the  entrance  angle  fg  ,  With  •  i 

4  (1-48) 

y 

Which  follows  from  (1,31?)  and  (1,0),  tho  bordorlino  of  stabi-  i 

lity  can  be  related  to  the  number  of  tooth  on  the  autter  z, 
so  that  the  ahart  is  valid  for  cutters  with  a  whatever  number 
of  tooth.  The  figure  describes  the  behaviour  for  tin  middle 

position  of  tho  table,  As  it  can  bo  seen,  tho  stability  ■ 

minima  in  the  stability  charts  In  Figure  1.28  o*-h  occur 
at  0°,  00°,  180°  and  2/0°  respectively,  At  thore  entrance 

angles  it  is  oquivnlcmt  for  tho  stability,  Whether  the  table  ■ 

1 a  moved  in  x-  or  y-direation,  booauso  for  n/D  "  0,B  the 
anglo  of  aontaot  is  just  90°. 

Considoring  tho  stability  charts  in  Figure  1.28  it  could  already  { 

bo  expected,  that  oonnldurablo  cUfforonooa  ouourod  in  the 
ohuttorfree  doapth  of  out  versus  tho  entrance  angle,  However 
in  practice  a  high  stability  in  only  important  and  uweful,  j 

if  it  can  be  reached  with  tho  installed  driving  power  and  j 

if  the  maximum  torsional  moment  must  not  be  excelled.  In  the  I 

program  RATTKH  2  those  limitation  can  he  calculated  too,  j 

When  calculating  and  plotting  only  the  minimum  stability  as  j 

performed  by  the  program  UATTKIi  2,  the  higher  stability  which  I 
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in  millinc?  operations  exists  in  Ronaval  between  the  unstable 
lobois,  is  not  taken  into  account  and  thus  cannot  turnod  to 
profit.  This  possible  increased  stability  is  preferably  intor- 
ostinfj  in  those  cases,  where  tho  machinin"  operation  -  for  ex¬ 
ample  at  a  workpiece  nlx'oady  mounted  -  must  be  pox’ formed  in  an 
unfavourable  configuration .  An  tho  evaluation  of  numerous  sta¬ 
bility  charts  for  faoo  milling  processes  showed  the  stability 
between  the  unstable  lobes  can  amount  up  to  about  200  %  of 


the  minimum  value,  if  the  ratio  of  tooth  mesh  frequency  f_  and 
ohattor  frequency  f^  satisfies  tho  condition 


ft  h  *  £ 
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Considering  the  ohattor  frequency  f.,  and  tho  value 


(t.40) 


which  iu'e  additionally  put  out  and  plotted  in  Viquro  t,  30,  the 


moat  unfavourable  rotational  upoedu  can  easily  bo  calculated 


by  meant)  of  the  equation  derived  from  oqu.  (1,.3d) 


6o  (1,1)0) 


In  praatiao  tho  increased  stability  can  be  turned  to  profit 
by  ahosiutf  rotational  upoodu  between  each  two  of  the  moot 
unfavourable  lipoods,  Am  it  udn  be  seen  in  tho  RiVoli  stability 
charts,  this  improvement  of  tho  outtlnn  performance)  is  specially 
possible  in  those  caset?,  where  only  ono  vibration  mode  id  pre¬ 
ferably  decisive  for  tho  chatter  behaviour,  If  the  frequency 
ratio  is  f a  /  i1^  <  0,535  or  if  more  than  ono  vibration  mode  in- 
fluonoo  tho  chatter  behaviour,  in  most  oases  tho  unstnblo  lobes 
overlap  in  such  a  manner,  that  tho  increased  stability  can  not  be 
ron lined. 


In  Vi  [turn  i,30  n  and  c  tho  minimum  orltinnl  doapth  of  nut  in 
plotted  unnimit  the  ontrivnco  ivni';le  for  D/lj  0,3  niul  }\/li  1  0,7, 
roiipoo lively ,  Per  D/D  .»  o ,  1,1  and  ifunurnlly  for  nil  11/1)^0,13  It 
in  possible  to  imuihlno  the  workpiece  with  table  travel  in  x~ 
as  well  an  in  y-dlrootion,  Vor  b/ll  "  0,7  and  puiuiruLly  ior  nil 
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1  j /D  >0,5  tho  curve  cuts  off  nt  a  certain  entrance  angle,  ainco 
the  Riven  brondth  of  the  workpioce  can  no  more  bo  machined  when 
exceeding  this  angle. 

Comparing  the  stability  charts  in  Figure  1.30  a,b,c  makes  ob¬ 
vious,  that  the  critical  deapth  of  cut  la  only  slightly  effec¬ 
ted  by  tho  D/1)  -  ratio  over  wide  ranges  of  the  entrance  angle 
for  instance  for  0°  90°  and  270°*  U  360°. 

This  result  means,  that  tho  utilization  of  the  installed  power 
beoomes  relatively  better  with  an  increasing  D/D-ratio  in  those 
rnngos,  whore  the  borderline  of  stability  lies  below  tho  limi¬ 
tation  duo  to  the  installed  power, 

Ill  the  shop  by  moans  of  such  representations  of  tho  stability 
behaviour  tho  critical  deapth  of  cut  to  bo  reached  can  easi¬ 
ly  be  predetermined  for  given  positions  in  the  working  range, 
given  tool  geometry  and  workpiece  material.  In  addition  to 
this  the  operation  personnel  can  get  a  good  qualitative  sur¬ 
vey  over  the  ohatterfroo  performance  of  the  machine.  Finally 
it  may  bo  stated  that  the  detailed  procedures  Will  onable 
the  maohihd  tool  users  to  improve  th©  utilization  of  the  ma- 
chines  considerably  in  many  oases, 
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2 ,  Chattel*  Specifications 


As  mentioned  above  the  purpose  of  the  research  work  performed 
on  the  e<r-otl  oo  "Chatter  Specifications"  wp.b  to  complete  the 
geometrical  acceptance  tests  (after  Schlesinger)  used  until 
now  since  over  30  years,  by  suitable  investigation  procedures 
and  specification  data  by  means  of  which  also  the  dynamic 
-  preferably  tho  chatter-  behaviour  of  machine  tools  oan  be 
tested  and  Judged  respectively.  The  necessity  of  dynamic 
acceptance  tests  of  machino  tools  results  from  the  customers 
demand  tor  minimum  cost  and  minimum  machining  time  at  pre¬ 
scribed  accuracy  and  surface  quality  of  the  workpieces.  The 
machine  tool  user  is  therefore  only  interested  in  a  procedure 
by  whioh  tho  cutting  performance  can  be  determined  as  perfect¬ 
ly  as  possible  under  consideration  of  tho  relative  displacements 
between  worltpioco  and  tool. 

At  the  Werkzeugmaaohinenlabor  in  Aachen  extensive  Investigations 
aimed  to  aoceptanae  tests  and  specif ioations  have  boon  carried 
out  especially  with  regard  to  the  ohattor  behaviour  jf  large  and 
heavy  maohino  tools.  Prior  conditions  for  this  wore  analytical 
investigations  ihto  the  chatter  fundamentals.  It  was  at  first 
neoessary  to  clarify  the  causes  for  l;he  obourvod  ohattor  pheno¬ 
mena  and  soaondly  to  work  out  methods  for  the  measurement  and 
the  calculation  respectively  of  tho  necessary  data. 

As  to  tho  question,  whiah  data  or  valuos  should  be  Ukiod  and  de¬ 
fined  for  specifying  ancl  judging  tho  chatter  behaviour,  when  ini¬ 
tiating  this  program  all  participating  research  groups  agreed,  that 
the  most  suitable  solution  was  to  tnko  tho  actual  critical  denpth 
or  width  of  out,  booauao  thin  is  the  value  only  being  interesting 
in  prautJou.  Thus  the  main  research  effort  was  at  first  directed 
at  the  problem  to  develop  the  nouesmary  techniques  for  determin¬ 
ing  these  values. 

Hereby  the  following  points  of  view  had  to  bo  const duvudi 


1,  A  minimum  number  of  values  snoulci  bo  sufficient  for  the 
judgement ; 

2,  The  measurement  procedure  and  the  calculation  method  should 
not  be  too  sophisticated  but  easy  to  understand; 

3,  The  required  effort  of  measuring  and  calculation  time  and 
the  tost-  and  evaluation  equipment  should  bo  as  minimum 
as  possible. 

In  tho  following  different  possibilities  resulting  from  tho  per¬ 
formed  research  work,  shall  bo  detailed  in  short  and  discussed 
with  regard  to  the  reliability  and  reproducibility  of  the  obtained 
result, s,  When  initiating  this  research  program  only  few  of  the  pa¬ 
rameters  having  influence  on  tho  chatter  behaviour  Wore  known. 
However  during  the  investigations  many  further  parameters  have  been 
found  out  of  partly  high  importance,  Tho  influencing  parameters 
for  one  part  already  quantitatively  for  tho  other  part  qualita¬ 
tively  known  until  now,  are  altogether  given  with  the  block  dia¬ 
gram  of  the  dynamic  manhilling  pi'ocosrJ  in  figure  2 , 1 ,  dome  of  thorn, 
especially  those  concernihg  tho  machine  dynamics,  au  for  in- 
atnnao  clamping  and  moving  of  tile  machine  elements  and  the  travel 
speed,  could  not  bo  investigated  before  the  development  of  tho 
measuring  procedure)  based  on  tho  random  noise  theory,  (see  [lo]  ). 
As  to  these  effects  only  few  maunuroiitolitu  oould  bo  performed  til 
now,  so  that  quantitative  relationships  being  generally  valid 
are  not  yet  availablo,  lSstamplud  for  the  inliuohuo  of  tho  chaugod 
position  of  tho  machine  elements  have  extensively  boon  disousnod 
in  [17,  l.lij  ,  invoutigationu  dotailad  In  [ill]  showed  tho  aonaidoru- 
bio  effect  of  tho  worltpibco  mounting  ,  Tho  dynamics  of  the  Work- 
pieuo  it so If  are  preferably  decisive  in  tho  higher  working  po¬ 
sitions  of  tho  iwohimiM.  Corresponding  moults  o;l’  investigations 
into  lurgu  portal  milling  mauliiueu  are  given  in  [loj, 

Mainly  throw  procedures  soum  to  be  possible ,  which  will  pro¬ 
vide  the  ncKuiwsnry  data  for  Judging  and  in  peel  fylng  the  chatter 
behaviour  of  machine  tools,  Two  of  them  provide  to  dot ermine  the 
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practically  interesting  cutting  condition*!,  those  prooeauros 
aro  Experimental  Chatter  Tests 

Theoretical  Chatter  Analysis 

The  third  possibility  ia  to  look  at  the 

Dynamic  Compliance  of  the  Machine 
at  the  cutting  point  because  there  ia  a  close  -  at  least 
qualitative  -  relationship  between  the  dynamic  compliances 
and  the  chatter  behaviour. 

2.1  Experimental  Chatter  Testa 

Considering  that  the  critical  deapth  or  width  of  out  is  mostly 
interesting  for  the  machine  user,  the  method  of  the  direct  me¬ 
asurement  of  the  stability  border  by  means  of  experimental  out- 
ting  tests  is  suggesting  itself.  The  main  advantage  of  maohi- 
ning  experiments  for  ohatter  analysis  of  machine  tools  is  the 
simple  performance  of  the  tests.  The  stability  border  can  be 
determined  by  altering  the  deapth  of  cut  and  observing  the  noise 
or  the  workpiece  surface,  More  comfortable  is  the  use  of  a 
vibration  pick-up  and  a  recorder  or  an  oscilloscope.  This  proce¬ 
dure  is  easy  to  understand  for  the  manufacturer  as  well  as  for 
the  user.  This  may  be  an  important  argument  under  oertain 
circumstances ,  if  trouble  raises  between  manufacturer  and  user, 
Further  important  advantages  are,  that  no  special  measuring 
device  and  qualified  personnel  are  required  for  performing  the 
tests. 

One  significant  disadvantage  is  t he  faot,  that  the  investigation 
requires  extremely  high  time  consumption,  if  numerous  combina¬ 
tions  of  the  influencing  parameters  must  be  taken  into  aooount , 
Tills  time  consumption  whereby  the  machine  must  bo  taken  nut.  of 
production  is  considerable  especially  at;  large  and  expensive 
machine  tools, 
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Additionally  several  problems  in  performing  tho  tests  must  be 
mentioned.  One  main  difficulty  of  chatter  experiments  1b  the 
fact,  that  the  stability  border  is  dependent  on  &n  external 
disturbance  force.  This  results  in  a  certain  hysteresis  of  the 
stability  border.  Going  from  the  stable  region  to  the  unstable 
yields  a  greater  critical  doapth  of  cut  than  coming  from  the 
unstable  region.  This  is  caused  by  tho  fact  that  near  the  Bta- 
bility  border  the  raise  timo  of  the  selfexcited  vibrations  de¬ 
pends  on  external  forces,  which  are  present  in  machining  n  Work¬ 
piece  surface  which  is  already  Wavy. 

Another  influence  factor  is  the  Workpiece  material,  In  (.  1 0 J 
results  of  experimental  chatter  tests  aro  given,  whore  the 
stability  border  is  plot  tod  versus  tho  diameter  of  a  turning 
WorkpAede,  Here  Considerable  differences  of  about  100  %  in  the 

critical  width  of  cut  wore  found  out  for  different  workpiece 
diameters.  Besides  to  thiB  the  influence  of  the  cutting  edge 
wear  must  be  observed  carefully,  HosUltB  of  experimental  chatter 
tests  oit  a  millinu  machine  £.  IB J  showed,  that  the  chatter  tests 
can  carried  out  only  in  a  certain  stage  of  tool  wear. 

finally  the  determination  of  the  critical  cutting  conditions 
by  means  of  experimental  chatter  teats  can  bo  regarded  as  the 
most  simple  procedure,  tho  reliability  of  which  however  is  li¬ 
mited  if  only  normal  time  consumption  shall  be  admitted.  The 
application  of  this  procedure  seems  to  bo  reasonable,  if  only 
small  spoctra  of  working  conditions  shall  bo  taken  into  account, 
for  Instance  at  special  or  single  purpose  machine  tools. 

Considering,  that  in  the  caHU  of  large  machine  tools  the  influ¬ 
ence  of  the  workpiece  upon  the  chatter  behavior  nan  not  at  all 
be  neglected,  the  difficulty  in  Hpeuifylng  the  above  muutiuiuid 
cutting  tustH  is  considerably  Increased.  Thus  It  can  bo  uiulur- 
stood,  that  thorn  are  many  obJocvlouM  against  the  general  appli¬ 
cation  of  cutting  lusts  for  specifying  the  nhattur  behavior  of 
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large  machine  tools,  though  the  procedure  itself  is  very  simple 
to  apply, 

,2,2  theoretical  Chatter  Analysis 

The  fundamentals  as  well  as  the  performance  of  tho  theoretical 
chatter  analysis  have  been  detailed  in  the  Quarterly  Technical 
Reports  Nr,  I,  II,  III  fli,  12,  13]  and  in  chapter  1  of  this 
report.  The  most  important  advantage  is,  that  by  means  of  the 
developed  computer  programs  the  stability  behavior  under  all 
interesting  conditions  oan  theoret ioally  be  determined  with  prac¬ 
tically  satisfactorily  acouracy  rationally  and  economically,  if 
the  characteristic  properties  of  the  machine  and  the  cutting 
process  are  known.  From  the  practical  point  of  view  It  is  ad¬ 
vantageous,  that  the  machine  must  only  be  available  for  the 
dynamic  investigation  whereas  the  evaluation  and  calculation  can 
be  carried  out  apart  from  tho  machine.  The  necessary  measuring 
time  at  the  machine  itself  will  furthermore  considerably  be  re¬ 
duced  by  applying  the  recently  developed  spectral  density  measu¬ 
rement  technique  developed  under  a  closely  related  contract  at 
the  Werksseugmaschinenlabor  der  Tl!  Aachen  too  [  laj  ,  By  means  of 
this  measuring  technique  it  Will  additionally  be  posuible  to 
take  the  machining  parameters  rotational  speed,  feed  rate,  moving 
of  elements,  clamping  a,s,o  into  account  during  the  dynamic  in¬ 
vestigation,  When  representing  the  chatter  behavior  in  the  form 
of  stability  charts  it  is  also  possible  to  calculate  and  to  plot 
tho  limitation  of  thu  cutting  performance  due  to  tho  installed 
power  and  duo  to  tho  maximum  torsional  moment,  hy  moans  of  these 
additional  values  it  is  easily  possible  to  compare  tho  chatter- 
free  performance  with  that  performance  which  the  machine  Is  de¬ 
signed  for,  and  therewith  gut  a  good  xml  clear  Judgemont  of  the 
practical  utilization  of  the  machine. 

Different  disadvantages  cannot  bo  let  out  of  cons  I dsral Ion ■  The 
first  one  In  due  to  tho  highly  sophist  inn  t  ml  nml  tints  not. 
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large  machine  tools,  though  the  procedure  itself  is  very  simple 
to  apply. 

2.2  Theoretical  Chatter  Analysis 

The  fundamentals  as  well  as  the  performance  of  tho  theoretical 
chatter  analysis  have  been  detailed  in  the  Quarterly  Technical 
Reports  Nr.  I,  II,  III  fll,  12,  13 ]  and  in  chapter  I  of  this 
report.  The  most  important  advantage  is,  that  by  means  of  the 
developed  computer  programs  the  stability  behavior  under  all 
interesting  conditions  can  theoretically  be  determined  with  prac¬ 
tically  satisfactorily  acouracy  rationally  and  economically,  if 
the  characteristic  properties  of  the  machine  and  the  cutting 
process  are  known.  From  the  practical  point  of  view  it  is  ad¬ 
vantageous,  that  the  machine  must  only  be  available  for  the 
dynamic  investigation  whereas  the  evaluation  ana  calculation  oan 
be  carried  out  apart  from  the  machine.  The  necessary  measuring 
time  at  the  machine  itself  will  furthermore  considerably  be  re¬ 
duced  by  applying  the  recently  developed  spectral  density  measu¬ 
rement  technique  developed  under  a  closely  related  contract  at 
the  Werkzeugmaschinenlabor  der  TH  Aachen  too  [id].  By  means  of 
this  meusuring  technique  it  will  additionally  be  possible  to 
take  the  machining  parameters  rotational  speed,  feed  rate,  moving 
of  elements,  clamping  a.s.o  into  account  during  the  dynamic  in¬ 
vestigation.  When  representing  the  chatter  behavior  in  the  form 
nf  stability  charts  it  is  also  possible  to  calculate  and  to  plot 
the  limitation  of  the  cutting  performance  due  to  the  installed 
power  and  due  to  the  maximum  torsional  moment.  By  means  of  these 
additional  values  it  is  easily  possible  to  compare  the  chatter- 
free  performance  with  that  performance  which  the  machine  is  de¬ 
signed  for,  and  therewith  get  a  good  and  clear  judgement  of  the 
practical  utilization  of  the  machine. 

Different  disadvantages  cunnot  be  let  out  of  consideration.  The 
first  one  is  due  to  the  highly  sophisticated  and  thus  not 
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generally  understandable  theory  the  whole  procedure  is  based 
upon.  Another  fact  is  the  complicated  measuring  and  evaluation 
procedure  as  well  as  the  very  expensive  test  equipment,  of 
which  the  operation  requires  still  qualified  personnel  up  to 
now.  Besides  to  this  other  difficulties  are  produced  by  the 
parameters  -  already  discussed  in  connection  with  the  cutting 
tests  -  influencing  the  chatter  border  during  machining. 

However  the  decisive  reason  which  is  qpposed  to  a  general  appli¬ 
cation  of  the  theoretical  chatter  analysis  is  the  fact,  that 
no  sufficiently  reliable  data  concerning  the  transfer  behavior 
of  the  cutting  process  are  available  at  this  time. 

Considering  the  advantages  and  disadvantages  the  above  detai¬ 
led  procedure  seems  to  be  the  optimum  possibility  for  future 
time,  if  the  chatter  behavior  of  machine  tools  shall  be  Judged 
by  specifying  the  actual  deapth  or  width  of  cut  for  a  large 
spectrum  of  machining  conditions. 


There  is  absolutely  no  doubt,  that  for  Judging  machine  tool 
chatter  behavior  the  actual  critical  width  or  deapth  of  cut 
are  the  most  suitable  and  interesting  values  from  the  practical 
point  of  view.  However  when  considering  all  the  mentioned  para¬ 
meters  (see  Figure  2,1),  there  raiseB  the  question,  whether 
it  is  reasonable  at  all  to  take  the  width  or  deapth  of  cut  as 
criterion  for  the  chatter  specifications  of  machine  tools. 

In  order  to  provide  a  test  as  simple  as  possible  it  seems  to 
be  advantageous,  to  let  the  parameters  of  the  cutting  process 
out  of  consideration.  In  this  case  no  longer  the  borderline 
of  stability  is  the  basis  for  the  Judgement,  but,  as  can  be 
seen  in  the  block-diagram,  the  dynamic  properties  of  the  system 
machine-workpiece.  If  the  compliance  of  the  workpiece  is  let 
out  of  consideration  too,  at  last  only  the  dynamic  properties 
of  the  machine  itself  must  be  regarded. 


2.3  Dynamic  Compliance  of  the  Machine 


The  fundamental  Idea  when  taking  only  the  dynamic  compliance 
of  the  machine  as  criterion  for  the  chatterfree  performance 
is  the  fact,  that  there  is  a  close-at  least  qualitative-relationship 
between  the  chatter  behaviour  and  the  mere  dynamic  compliance  of 
the  machine,  To  test  these  properties  is  relatively  simple,  but 
it  must  be  verified,  what  judgement  about  the  chatter  behavior 
can  be  made  by  considering  only  the  dynamics  of  the  machino. 

As  it  is  simple  to  understand  this  test  provides  not  at  all  an 
absolute  judgement  of  the  machine  tool.  However  as  to  the  largo 
machine  tools  it  seems  to  bo  highly  questionable,  whether  such 
a  test  is  still  reasonable,  after  the  invest igations  performed 
till  now,  have  shown,  what  decisive  influence  upon  the  chatter 
behavior  is  often  caused  by  the  compliance  of  the  workpiece 
(see  [  10,  15 J  ).  To  this  it  must  bo  replied,  that  first  of 
all  the  machine  itself  must  fulfill  certain  prior  conditions 
concerning  the  dynamic  stiffness  in  the  critical  vibration 
modes . 

Considering  for  instance  a  portal  milling  machine  the  most  cri¬ 
tical  vibration  modes  are  the  bending  modes  of  the  portal. 

Presuming  now  the  workpiece  to  be  rigid,  it  can  be  stated,  that 
the  chatter  behavior  of  a  machine  A,  which  in  the  two  vibra¬ 
tion  directions  of  the  portal  has  by  far  higher  stiffness  as 
a  machine  B,  will  be  better  than  the  chatter  behavior  of 
machine  B.  By  this  simple  example  it  shall  only  be  demonstrated, 
that  by  means  of  the  dynamic  properties  of  a  machine  a  relative 
judgement  between  machines  of  the  same  typo  is  possible  at  all. 

On  the  other  hand  it  must  be  noticed,  that  considerable  uncer¬ 
tainty  may  be  exist,  when  applying  this  method.  It  will  not 
at  all  be  possible,  to  find  out  differences  of  20  to  30  %  bet¬ 
ween  machines,  however  it  should  be  possible  to  ascertain  for 
example  a  machine,  which  is  about  300  %  worse  than  other  ones. 
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Thus  it  can  be  said,  that  this  method  will  provide  a  good  en¬ 
gineering  judgement  about  the  dynamic  behavior  and  therewith 
about  the  chatterfree  performance.  However  before  generally 
applying  this  method  extensive  dynamic  and  chatter  testB  will 
still  be  necessary  into  the  problems,  what  vibration  modeB  are 
decisive  for  the  chatter  behavior  at  the  different  types  of 
machine  tools  and  in  howfar  the  dynamic  and  the  chatter  be¬ 
havior  vary  due  to  the  changing  machining  parameters  as  for 
instanco  the  position  of  the  elements  within  the  working  range, 
the  movement  of  the  elements  and  the  clamping  conditions. 


Summarizing  the  considerations  concerning  the  chatter  specifi¬ 
cations  it  can  be  stated,  that  now  after  having  systematically 
investigated  and  clarified  the  fundamental  relationships  being 
valid  in  machine  tool  chatter  three  procedures  to  test  and  spe¬ 
cify  respectively  the  chatter  behavior  of  machine  tools  aro  now 
available  on  principle.  However  the  above  discussion  of  these 
procedures  showed,  that  further  practical  cutting  tests  as  well 
as  dynamic  investigations  are  necessary  before  one  of  these 
procedures  can  be  proposed  to  be  the  most  suitable  one  in  gene¬ 
ral  and  before  general  purpose  maohine  tool  chatter  specifications 
can  be  established  respectively. 


3.  Structure  Analysis 


3.1  Digital  Computation  of  Machine  Tool  Elements 
andStru  cture  s 

For  the  improvement  of  the  dynamic  behaviour  of  machine  tools 
it  is  an  important  point  oi  view  to  develop  analytical  methods 
for  the  calculation  of  the  dynamic  properties  already  in  the 
status  of  design.  As  the  mathematical  treatment  of  a  complete 
machine  tools  system  is  an  extremely  complex  problem,  due  to  the 
influence  of  joints  and  the  interdependence  of  the  single  ele¬ 
ments,  thoso  investigations  must  start  from  a  separate  analysis 
of  the  different  components.  The  next  stop  will  then  bo  the  connec¬ 
tion  of  these  elemonts. 

All  the  analytical  methods  being  developed  have  been  programmed 
and  tested  by  comparing  measured  an  computed  results.  The  methods 
are  described  in  u  form  which  makes 'it  possible  to  reproduce 
them,  easily  on  any  computer. 

3.1.1  Computation  of  Spindle  Systems 

Sevoral  investigations  of  heavy  machine  tools  have  shown  that 
the  characteristics  of  the  working  spindle  are  of  great  impor¬ 
tance  for  the  working  accuracy  and  the  chatter  behaviour.  It 
was  stated  that  under  certain  conditions  up  to  one  third  of  the 
total  deflection  is  given  due  to  spindle  bending.  (Fig.  3.1) 

The  computation  techniques  for  spindle  systems  known  from  the 
literature  so  far  usually  require  considerable  simplifications 
concerning  the  geometrical  shape  or  the  boundary  conditions  e.g. 
the  arrangement  of  the  bearings.  In  the  following  a  method  will 
be  presented,  which  is  applicable  to  spindles  with  any  arrangement 
and  number  of  bearings. 

Usually  the  cross-sections  of  a  spindle  and  with  this  the  moment 
of  inertia  I  (x)  and  the  mass  m(x)  will  vary  along  the  spindle- 
axis.  According  to  Timoshenko  the  behaviour  of  a  spindle  element 
ns  shown  in  Figure  3.2  is  described  by  tho  partial  differential 
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aquation 


(3.1) 


A  complete  solution  of  equation  (3.1)  is  generally  impossible. 
But  under  certain  conditions  a  sufficiently  aorrect  approxima¬ 
tion  can  be  obtained  when  using  electronic  computers.  For  the 
solution  of  the  problem  by  means  of  a  digital  computer  the 
temporal  deviation  is  made  equal  to  zero.  Then  equation  (3.1)  is 
transformed  into  an  ordinary  differential  equation 
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This  can  be  solved  for  an  elastic  beam  with  many  different  cross- 
sections  and  several  flexible  supports  by  a  method  using  matrix 
algebra.  The  following  results  can  bo  obtained. 


a)  static  behaviour  of  the  system  spindle-bearing, 

b)  natural  frequencies  of  the  systems. 

For  the  computation  the  spindle  is  divided  into  separate  sec¬ 
tions  which  may  have  different  lenghts.  In  order  to  achieve 
a  high  accuracy  as  many  sections  as  possible  should  bo  chosen. 
The  state  of  forces  and  deformation  7^:  is  described 

by  the  vector 

y 
r 

M 
Q 

The  general  solution  of  equation  (3,2)  for  a  ny  section  "k"  of 
the  spindle  (see  Fig.  3.2)  is  obtained  by  four  fold  integra¬ 
tion  for  the  region  0  *  xk  *  l[t  [ill. 

After  that  the  relationship  between  any  two  interesting  points 
of  the  upindle  may  bo  expressed  in  matrix-notation. 
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Additional  influences,  such  as  flexible  bearings,  uro  takon 
into  consideration  by  a  so-called  point-matrix  p^.  This  opera¬ 
tion  can  be  explained  by  means  of  Figure  3.3.  The  status 
dlreotly  on  the  loft  side  of  the  point  b,  which  is  cullod  b^ 
starting  from  tho  status  on  the  right  side  is  obtained  by 
the  calculation 


tyhi  ■  Xk  1  ^  m 


(3.4) 


respectively,  where  is  the  so-callod  transfer  matrix. 


External  forces  (j  or  moment  M  and  alterations  of  the  diameter, 
which  cause  an  alteration  of  the  load  q,  are  taken  into  account 
by  defining  the  vector  in  a  more  general  form 


h 


»• 


>'« 

r* 

*%*  * 
Q*  t  a 

vs  •'  vL 


(3.5) 


The  values  and  f ^  of  point  1,  which  are  unknown  so  far, 
are  obtained  from  a  linear  system  of  two  eauations  in  depen¬ 
dence  on  tho  coundary  conditions,  M  “  0  and  Q  =  0  [  1  l/. 


If  the  initial  status  vector  ^  is  known  the  vector  for  a 
given  point  can  be  determined  in  the  following  way, 

t. 

-£s  -  fjX  -*  /Jt 

:  (3.6‘) 
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In  Figure  3.4  the  measured  and  the  calculated  bending  lines 
of  a  lathe-spindle  loaded  at  the  spindel-nose  by  a  static 
force  of  2200  lb  (1000  kp)  aro  compared.  The  correspondence 
between  the  measured  and  calculated  deflections  is  quite 
good. 

In  order  to  determine  the  natural  frequencies  a  mass-matrix, 
similar  to  the  point-matrix,  is  introduced  into  the  program 
for  the  static  bending  line  [l3j.  Assuming  an  arbitrary  deflec¬ 
tion  of  the  spinclle-noso  the  solution  is  found  in  an  iterative 
computation  process.  In  the  natural  modes  both  the  force  0 
and  the  moment  M  at  the  end  of  the  spindle  must  be  equal  to 
zero.  The  frequencies  where  this  condition  is  fulfilled  are 
calculated. 

By  aid  of  this  method  given  in  form  of  a  block  diagram  in  Fig.  3.5 
the  natural  frequencies  and  according  mode  shapes  may  be  found. 

If  there  are  only  the  natural  frequencies  of  interest  the  effort 
can  be  decreased  rapidly.  [1.3 j 

Figure  3.6  shows  the  resonance  frequency  of  a  lathe-spindle 
with  a  workpiece  clamped  in  the  three- jaw-chuck  as  a  function 
of  the  front  bearing  stiffness.  It  can  be  seen  that  tho  dyna¬ 
mic  behaviour  of  the  system  is  described  quite  exactly  by  the 
chosen  method. 

Furthermore  in  a  test  series  the  natural  frequencies  of  five 
machine  tools  of  the  same  type  and  the  same  type  of  spindle 
systems  have  been  measured.  A  remarkable  scattering  of  tho  na¬ 
tural  frequencies  has  to  be  stated  which  cannot  bo  explained 
completely.  The  age  differences  of  the  single  machines  do  not 
show  a  definite  influence.  Using  tho  digital  computer  program 
and  the  mathematical  model  of  the  spindle-bearing  system  shown 
in  Fig.  3.7  the  natural  frequencies  have  been  calculated.  A 
comparison  between  t!  eorefical  and  experimental  results  -  see 
Fig.  3.8  -  yields  deviations  between  0  and  -16,5  %.  On  tho 
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average  r  difference  of  -  7,25  %  is  found.  If  this  is  com¬ 
pared  to  the  scattering  of  the  actual  natural  frequencies, 

+  11,1  %,  the  usefulness  of  tho  computation  method,  even 
if  so  complex  structures  as  the  working  spindle  of  a  boring 
mill  are  concerned,  becomes  evident. 

When  applying  the  chosen  procedure  of  computation  the  ncouracy 
of  the  results  is  the  higher  the  better  the  mass  distribution 
is  approximated  within  tho  mathematical  model.  However,  from 
this  fact  it  cannot  be  concluded,  that  the  accuracy  could  bo 
improved  at  any  rate  by  increasing  the  number  of  mass  elements 
as  at  the  same  time  the  effect  of  rounding-errors  become  more 
significant.  Duo  to  these  opposed  influences  an  optimum  for  the 
number  of  mass  elements  will  exist.  In  the  present  aa.se  the 
computer  program  is  suitable  for  a  separation  into  28  mass 
elements  and  no  detrimental  effect  due  to  rounding  errors 
was  observed. 

3.1,2 A  numerical  method  for  calculating  gross  section  values 

In  this  chapter  a  numerical  method  will  bo  described  to  cal¬ 
culate  the  cross  section  values  of  an  arbitrary  given  poly¬ 
gonal  figure.  By  help  of  this  method  the  following  values 
can  be  calculated: 

Shear  Area, 

Neutral  Axis, 

Area  Momonts  about  given  axis, 

Inertia  Moments  about  given  axis, 

Inertia  Moments  about  Neutral  axis, 

Counter  clockwise  angle  of  rotation  of  principle  axis, 

Inertia  Moments  about  principle  axis, 

Polar  Mass  Moment  of  Inertia, 

Transverse  Mass  Momonts, 

Weight  of  element. 
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All  these  values  are  needed  as  inputdata  for  the  programs 
to  calculate  the  static  and  dynamic  behaviour  of  columns 
and  frames  of  machine  tools  as  described  in  chaptor  3,1.3, 

The  calculation  of  shear  area  values  of  single  or  inerrupted 
areas  as  shown  in  Fig.  3. Dean  be  done  by  the  Gauss  Integral 
Law  for  areas,  called  Groen's-Law  as  well. 

In  the  most  general  form  this  law  may  be  written 

R  •  Jf  f(*,v)  clxcly  f  F(x,y)  ,  (3.7) 

3  c. 

On  tho  other  hand  for  a  given  area  in  the  x-y-plane  the  shear 
area  value  can  be  calculated  by  the  following  formula  too: 

R  -  ff  ctx  d  y  (3.C) 

From  this  formula  it  can  easily  be  obtained  that  F  (x,y)  in 
formula  3,7  has  to  bo  equal  to  1  in  this  case.  Thus  F  (x,y) 
in  the  right  hand  part  has  to  be  equal  to  x: 

F.  ,  i 

dx 

Now  formula  (3.7)  may  also  be  written 

F  “  Jfd*  c(y  *  j>xdy  (3,10) 

Taking  into  account  that 

$  X  c/y  ‘  ~  f  y  dx 

c  c 

it  can  be  written 

*  %  f(xdy-ydx)  (3.11) 

c 

This  equation  can  be  mplitted  In  general.  Pro.jumi.ng  all  the 

areas  to  be  calculated  being  given  as  polygonal  figuros  this 
formula  can  be  subdevided  into  a  lot  of  segments  where  each  of 
them  is  sot  up  along  the  straight  line  be two on  any  two  given 
polygon  points.  Thus  for  equation  3,11  tho  following  numerical 
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formula  may  be  written 
R*  %  ± 

L-  t 


l*i  (\*. 


X  )  ~  K  (■*!'■> "  X  jJ 


or  in  a  more  compact  form 

n 

h  •  Y*  2  Lti  -  y{  i,w  J 


(3.12) 


(3.12a) 


In  this  formula  the  index  i  is  following:  the  curve  C  in  mathe¬ 
matical  positive  way  as  shown  in  figuro  3.10  and  for  i  i>  n  point 
(n  i-l)  is  equal  to  the  starting  point. 


Tho  formula  as  given  in  2.12  a  can  easily  be  programmed  to  calcu¬ 
late  tho  value  of  the  cross  section  area. 

The  next  step  is  to  calculate  the  areamoments  and  tho  noutral 
axis. 

Therefore  always  two  neighbourpoints  with  the  subscripts  i  and 
i  +  1  and  the  starting  point  of  tho  system  of  coordinates  are 
assumed  to  be  the  three  corner  pointsof  a  triangle.  For  each  of 
the  n  triangles  which  are  possible  the  neutral  y-axis  for  example 
oan  bo  given  as 

y  Yi.*+Yt' 

'i  -  3 

The  area  of  each  of  them  can  bo  calculated  as 

Ct  x  7*  (  *4  Y -Vv  “  X  x  lf-4  ) 

Thus  the  area  moment  of  such  a  trianglo  becomes 

5r  -  “  -y*’  %.  ('.  -K  )  •  %  ■  ( '  Y' ) 

As  long  as  it  is  possible  to  surround  those  triangles  in  a  mathe¬ 
matical  positive  way  with  respect  to  the  ascending  numberod  points 
the  values  for  s  become  positive  otherwise  negative  (Fig.  3,10). 
The  sum  of  all  of  thorn  is  equal  to  the  area  moment  of  the  whole 
cross  section.  About  tho  x-axis  it  can  be  written 

S.  '  Vi  f-  L(x  ■  X  )U  4  )  ] 

1-4  J 
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(3.13) 


and  about  the  y-axis  in  the  sane  way 

5*  %,  %  C(ft  ft**  -ti.iYiXti'-*.-**)  2 

'  4  ■«# 


(3.14) 


From  these  equations  now  the  distances  between  the  given  x-y 
axis  and  the  neutral  ones  can  be  calculated: 


n  ' 


s, 

0 


(3.15) 


Quite  similar  to  the  results  as  shown  until  now  the  formulas 
for  the  moments  of  inertia  about  given  x-  and  y-axia  can  be 
found 

n 

\  •  7*  -  ft**  Vi  )-l(y,  *  y;+.<)X-  Yi-Vii.4  ]]  (3.1G) 

l 

and 

1ym  %3  X  L  (  ft  ft**  ~  ft'**  ft  )  L  '  X;  '  ii+4  J  J  (3,17) 

1*4 


For  tho  moment  of  deviation  I„„  about  the  given  axis  it  can  be 

xy 

written 


0xy‘  Zil '{(XiK-'4-Xii<Yi}[(*i**«*)(Yi*yi'-4)T  X  (ft  j 


(3.18) 


All  these  formulas  are  easily  to  bo  programmed  and  the  further 
existing  characteristics  of  any  given  polygonal  oro3S  section 
can  bo  calculated  on  the  basis  of  those  values. 


Tho  moments  of  inortia  being  calculated  about  the  given  coordi¬ 
nates  can  be  expressed  about  the  neutral  axes  by  help  of  Steiner 
Law, 
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(3.19) 


s 


Using  these  expressions  the  angle  of  rotation  between  the  noutral 
axis  system  and  the  system  of  principle  axes  (Fig,  3.11)  can 
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be  calculated 


tnn(Jty)-  (3.20) 

J  w  -  Jr  <*■ 

Furthermore  the  principle  moments  of  inertia  about  the  ?-  and 
^  -axis  are  given  by  the  following  formulas 

If  «  Iff.rJxt  ,  j-  fa. 7  JjSt  cotUlr)  -  3*,s-  sinf2  V)  I 

Op  .  -  JT  caste*)  -  •  tlnd'S)  j  (3.21) 

On  tho  basis  of  the  cross  section  values  as  described  above  it 
is  easy  to  obtain  tho  values  for  tho  weight  of  an  element  of 
the  lenght  L  (Fig.  3.12)  and  the  polar-  and  transverse  mass 
moments  of  inertia. 

In  Fig.  3.13  a  block  diagram  of  the  program  is  given  for  cal¬ 
culating  all  the  values  mentioned  above.  The  applicability  of 
tho  method  is  shown  for  the  cross  section  of  a  machine  tool 
column  as  given  in  Fig.  3.14.  The  results  of  calculation  as 
they  are  put  out  from  the  computer  are  shown  on  Fig.  3.15. 

3. 1.3  Calculation  of  the  Static  and  Dynamic  Behaviour  of 
Machine  Tool  Elements  and  Frames 

110  Fl.g«  3.1  the  f lux-of-force-analysis  for  a  heavy  boring-  and 
milling-machine  is  shown.  The  result  of  this  analysis  is  typi¬ 
cal  for  most  of  these  machines.  It  can  be  seen  that  about  1/3 
of  the  total  deflection  results  form  the  spindle-system. 

Methods  for  calculating  the  static  and  dynamic  behaviour  of 
this  part  have  been  presented  and  applied  to  some  spindle 
systems.  As  the  next  stop  a  method  for  calculating  the  static 
deflection  of  a  machine  structure  -  reduced  to  a  throo-dimen- 
sional  beam-system  -  will  bo  described. 

3 . 1 . 3 . 1  Calculation  of  Static  Deflections 

'J'ho  method  described  here  is  basod  upon  tho  possibility  to  re- 
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duce  a  machine-tool-structure  to  a  three-dimensional  system 
of  uniform  beams.  This  is  possible  in  the  very  mcnt  cases.  An¬ 
other  limitation  is  Hooke's  Law. 

After  the  maohine-tool-structure  has  been  reduced,  the  diffe¬ 
rent  elements  have  to  be  named.  The  denomination  can  be  seen 
in  Fig.  3.16.  All  the  elements  are  marked  by  the  letters  of 
the  alphabet  beginning  with  "a". 

Internally  the  beginning  of  eaoh  element  is  called  "1",  the 
end  "2",  The  ends  "2"  of  eaoh  element,  which  are  the  conneo- 
ting-Joints  between  two  elements  at  the  same  time,  are  indica¬ 
ted  with  the  capital-letters  of  the  element's  names.  The  de¬ 
flection  of  a  point  A,  for  instance  in  the  x-direotlon  due  to 
a  load  applied  at  point  0,  will  be  called  x^  q.  The  rotations 
are  indicated  in  the  same  way. 

In  tho  following  two  systems  of  coordinates  are  used.  One  is 
a  right  handed  general-system  for  coordinating  all  the  members 
the  axis  are  x',  y'  and  z',  On  the  other  hand  each  of  the  ele¬ 
ments  has  its  own  system,  the  member-system,  which  is  a  right 
handed  one  too.  The  axis  of  the  member-system  of  an  element  i 
are  called  x^,  yi,  and  z ^  where  x^  is  always  the  lenght-axis. 

The  sign-declaration  (Fig.  3.17)  is  the  one  of  a  right  handed 
system  of  coordinates  for  the  deformations  as  well  as  for  the 
loadings. 

Each  of  the  elements  is  regarded  to  be  rigidly  supported  at 
point  "1".  The  real  boundary-conditions  will  be  introduced 
later.  For  such  an  clement,  loaded  at  "2"  the  defleotions  and 
distorsions  may  be  expressed  with  respect  to  Hookes 's  -  Law  as 
•f,  L/ftr  o  0  0  a  o  f>t 

Jy  0  L>/uj.  <?  o  a  LXsjt  t°, 

*  -  a  O  LJAsjy  0  -  0  I  I  I 

0  0  0  l/fir  0p  (3.22) 

//  0  0  0  Vejy  0 

0  Ly*£Ji  0  0  0  VfJt  m. 
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or  d  •  F.  f 

(3.22a) 

The  F-matrlx  is  symmetrical  (Maxwell-Betty  reciprocal  theorom). 

Alter  the  relationship  between  force  and  deflection  has  been 
demonstrated  for  one  arbitrary  element  of  a  system,  the  tran¬ 
sition  to  a  complete  system  will  be  explained.  The  main-prob¬ 
lem  in  a  complex  system  in  space  is  the  so-called  shifting  of 
forces  and  moments  from  the  acting  point  into  all  the  other  ones 
for, computation. 

By  looking  at  the  compatibility-conditions  between  the  ends  "1" 
and  "2"  of  an  arbitrary  placed  and  formed  element  (Fig.  3.18) 
therefore  the  following  relationship  can  be  obtained. 


(3.23) 


By  this  equation  a  force  can  be  shifted  in  the  moat  general 
case  within  a  system  and  the  ends  "1"  and  "2"  may  be  any 
two  point  I  and  S: 

r%  *  •  f*,  .  <r  ( 3.23  a  )  ;  ^  s  •  f*i  •  <0-  (3,23b) 
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The  total  deflection  as  well  as  the  torsion  of  any  system 
joint  can  be  subdivided  into  two  components! 


(3.24) 


a)  the  part  due  to  the  elastic  deformation  of  tho  rogarded  olo- 
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ment  itself,  caused  by  a  load  acting  anywhere,  called  member 
deformation  ©  and 

b)  the  part  due  to  the  deformation  of  the  other  elements  (in 
accordance  to  the  general  denomination,  which  is  called  rigid 
body  displacement  H*. 


Thus  for  the  total  deflection  and  troaion  it  may  be  written 


•r/ 

of y 

ey 

Jz 

<?* 

** 

A 

. 

r-\ 

•  i 

or 


3.  -  £ 


For  the  member  deformation  may  bo  written 


3* 


(3.25) 


£  .  F  •  yo  (3.26) 

To  find  3*  we  have  to  consider  a  completely  rigid  member  with¬ 
out  any  elastic  deformation.  Its  end-deformations  shall  be  3^ 
and  32  due  to  the  forces  p^  and  p2.  For  the  total  virtual 
work  being  done  by  tae  applied  forces  it  may  bo  wrtitten 


■r4 1  .  <4  -f/'  <4  -  ^ 

(3.27) 

and  substituting  from  equation  (3.4) 

-  rS  ■  h*-  +  r 3 '  cTi  -  &■ 

(3.27a) 

it  follows 

(3.2U) 

With  equation  3.2H  any  existing  deformation  at  point  lean  be 
related  to  point  2  and  particularly  the  rigid  body  deformation 


<7/  - 


(3 . 2Ua ) 


Substituting  a  from  oqu,  3.26  and  3*  from  3.2Hu  the  total  Uo- 
i’ormation  can  he  expresses  for  any  given  element  k  between 


-  !i7  - 


the  joints  I  and  K: 


'•C  ~  Fk 


n  ,  ///.  <3-29> 

In  the  previous  chapters  only  systems  consisting  of  elements 
without  any  relativ  angles  between  them  were  discussed.  To 
handle  systems  where  arbitrary  angles  between  the  members  may 
occur  two  possibilities  exist: 


a)  The  forces  and  deformations  can  be  transformed  from  one 
member-coordinato-system  into  another  one  at  each  joint  or 

b)  all  the  forces  and  deformations  of  any  element  are  trans¬ 
formed  into  a  general  coordinate-system. 


The  last  possibility  is  chosen  for  the  following  investigations. 
For  all  the  acting  foroes  and  moments  the  general  transorma- 
tion  matrix  T  can  be  obtained: 


r* 


csjYOJ  cas(x',yk)  cv,  fr'ti.)  0 

cojfr.'i co s  (*]  yk  )  cc>i(y!tk)  p 

cpsU'hJ  cos(z',Yk)  co%(t\zk-)  o 


o 

0 

0 


o 


0 

o 


(3.30) 


0  0  cos(t;*k)  eos(/,yu) 

0  0  e»s(y'/K)  cos  (y'y¥)  <°i(y',tk) 


0 


0 


0  COS(x‘/tk)  SOS  COtft't") 


The  transformation  matrix  T  is  orthogonal.  Thus  it  may  bo 

written : 

jn  <*  -4 

7  "  r  (3.31) 

The  relationship  between  t,-.e  deformations  in  the  member  system 
and  the  general  system  can  bo  found  in  the  same  way  as  shown 
for  tho  loadings 

r~'  “  T ■  *  (3.32) 
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Buoatituting  d  by  means  of  equations  (3.22a)  it  can  be  written 


<x  •  r  -  r  ■  i S' 

The  H-matrix,  which  allows  to  shift  loadings  within  the  system 
can  be  transformed  in  tho  same  way; 

H  '  -  r-  h  ■  r'i 

A  vector  p  will  bo  introduced  which  yields  all  tho  forces  at 
the  ends  "2"  of  the  elements  in  member-coordinates 


A* 

nc 

\rfllh 

%  * 

In  the  same  manner  a  vector  cm  is  defined  which  contains  all 
the  member  deformations  at  the  ends  "2"  of  the  whole  system 


Now  the  relationship  between  oju  and  p  due  to  the  elasticity 
of  the  elements  can  bo  written  .In  accordance  with  equation  (3.2G) 

o  c  o  ■ 

0  O  •  « 

‘  -  (3.37) 

■>  o  o  1 


Before  the  rigid  body  deformation  or  boundary  conditions  will 
be  included  the  formulas  obtained  so  far  will  be  transformed 
into  the  general  coordinate  syslc;.  . 


It  may  be  written: 

/-i  -  T-  /%  -  7"  //  • 


(3.19) 
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The  equilibrium  conditions  of  the  joints  in  Fie.  3.19  can  now 
be  formulated  in  general  coordinates: 


* 

-  ^ 

’  r> 

"  r*  •  f\€ 

(3.28) 

or  in  matrix  notation 

ra  -  T*  ■  h\ 

*  r* 

0  P 


P 
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K  ■  z-4 

rc 

i 

i 

(3.28a) 


or 


c  '  /v 


(3.28b) 


The  matrix  C  Is  called  connection  matrix. 


In  addition  the  following  relationship  between  the  member-  and 
total-deformation  can  be  obtained 

m  C*  •  J)'  (3.29) 

Substituting  in  equations  3,28b  by  the  formulas  3.26  and 
3.29  it  may  be  written 

p'  -  (3,30) 

or 

r  <?"*■**  ■€'*]■  P'  (3.31) 

This  equation  represents  the  relationship  between  loadings  and 
deformations  of  an  arbitrary  system  in  the  most  general  case. 

In  equation  3.31  the  11-matrix  occur  in  momber  coordinates. 
Usually  it  is  oasier  to  take  the  length-olements  for  the 
11-matrix  directly  out  of  a  drawing  in  goneral  coordinates. 

In  order  to  introduce  theso  values  into  the  C-matrix  the 
11-matrices  have  to  be  substituted  by  the  H-ma trices. 
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Another  simplification  ife  to  indicate  the  H-matrlces  by  means 
of  the  joint  names  as  for  instance 

Furthermore  H^-matrlces  aro  introduced,  which  are  equal  to 
the  unit  matrix,  in  order  to  get  more  systematology. 


Now  for  the  final  matrix,  given  in  equation  3.31,  it  can  be 
written 

' - - - - 7  I 
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Substituting  each  element  of  the  final-matrix  (oqu.  3.32) 
by  a  from  like  for  instance 


Hp,  T.  Fa  ■  Tu*  ■  h/M 

equation  3.24  can  be  simplified 
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(3.33) 


The  applicability  of  this  method  can  be  obtained  from  the  com¬ 
parison  of  the  moasured  and  computed  deformation  of  a  compli¬ 
cated  machine  tool  column  as  given  in  Fig.  3.20,  The  computed 
results  are  shown  in  Fig.  3,21  as  they  were  pointed  out  by  the 
computer.  For  point  6  for  a  computed  deflection  of  7,13  ^um 
(0.282  thou)  obtainos.  This  results  corresponds  very  good  with 
the  measured  value  of  7.25  yum  (0.2B6  thou). 

The  above  mentioned  method  can  be  used  for  calculating  multiple 
pinnod  systoms  too.  Therefore  the  load  vector  of  equation 

d-[F\  P 

has  to  includo  the  reaotion  loads  at  the  additional  pinned 
points.  These  loads  can  easily  bo  found  by  the  principle  of 
superpositioning  ^16  j  . 


G2 


3. 1.3, 2  The  Calculation  of  olgonfroquenc i o s 


The  differential  equation  for  an  undamped  one  mass  system  Is 
/Yl  K  +  k-x  “  O 

In  the  same  way  it  is  possible  to  describe  the  dynamio  beha¬ 
viour  of  a  system  of  soveral  masses  connected  by  springs: 

[/vjy  +•[/<']')<'  =  O  (3.34) 

In  this  equation  the  expression  J.K]  roprosonts  the  stiffness 
of  the  whole  system  and  [M is  tiie  so  called  mass  matrix. 

The  inverse  of  [ K  J  is  called  flexibility  matrix  of  the  system. 
The  development  of  this  matrix  has  already  been  deserxbod  in 
the  Quarterly  lleport  Nr.  4  for  single  systems  and  in  Quarterly 
Report  Nr,  6  for  multiple  pinned  ones  [14,16], 


The  mass  matrix  in  general  is  a  diagonal  matrix  of  the  form 


in  A 

,mL  O' 

rrl  I 


0 

y>h, 


(3.35) 


for  a  given  system  consisting  of  n  elements  and  n  lumped  masses. 
As  for  any  element  in  space  six  degrees  of  freedom  are  taken 
into  account  each  of  the  elements  of  tho  mass  matrix  consists 
of  a  diagonal-matrix  itself. 


To  find  out  the  oigenfroquencies  and  modoshapes  of  a  given  sy¬ 
stem  described  in  short  hand  notation  by  equation  3,34  it  can 
bo  assumed 


x.. 


(3.4) 


Introducing  this  expression  into  equation  3,1  it  can  be  written 


\tl}  lot  X,  *  [k]  X. 


o 
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and  with 


Ik)  -  [f\ 

[f!  M  -  ftf] 

and 

i^r 

this  formula  becomes  the  most  general  form  of  an  eigenvalues- 
problem 

([fl]  -  l-  l:  )  '  X  '  O  <3*35) 


After  the  n  oigcnfre  quencios  of  a  systom  consisting  of  n  ele¬ 
ments  and  lumped  masses  have  boon  found  the  next  step  Is  to 
determine  tho  according  eigenvectors  Y,  . 

Therefore  the  system  of  equations  3,35  has  to  be  solved  for  all 
the  eigenvalues : Ar  : 


l\fj]  -A,£  -  \B),  X-  -O 


(3.3(3) 


After  the  eigenfrequonoies  and  modeahapes  liavo  boon  calculated 
it  is  possible  to  calculate  tho  response  locus  for  a  chosen 
degree  of  freedom  [  B  J  ,  tho  damping  values  which  are  noedod 
therefore  have  to  bo  assumed  previously. 


To  find  out  the  correlation  between  computed  and  measured  re¬ 
sults  the  machine  tool  colomn  of  Fig.  3.32  was  oxcitod  in  the 
x-  and  y-direction.  The  test  rig  is  shown  in  Fig.  3.23.  The 
measured  and  computed  first  eigenfrequencies  for  both  direc¬ 
tions  are 


direction 

measured 

1  frq. 

computed 

1  frq. 

error 

X 

326  Hz 

335  Hz 

+  2,76  % 

y 

300  Hz 

323  Hz 

+  7,66  % 

Finally  the  frequency  response  loci  for  both  directions  were  cal¬ 
culated  as  already  described  before  on  the  base  of  the  so  far 
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computed  eigenfrequenci.es  an  according  modes  shapes. 

The  results  for  both  directions  are  shown  in  fig.  3.24  und  3.25 
as  they  were  printed  out  from  the  computer  for  exciting  and 
picking  up  at  the  top  of  the  column.  The  correlation  between 
the  measured  and  aomputed  ampJLtudo  for  the  above  mentioned 
first  frequencies  can  be  obtained  from  the  following  table: 


aplitude  for  the 

first  bending  eigenfrequoncy 

direction 

measured 

computed 

error 

x 

3,4  ^um 

3,01 

-11,5  % 

y 

3,3  /Um 

3,33  /am 

0  % 
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3 , 2  Thoorotical  Analysis  of  Loading  and  Doformatlon  of 
Columns 


For  analysing  tile  elastic  dcformntion  oi’  columns  some  simpli- 
fioations  have  to  bo  assumed  both  for  loading  and  sti'osa  distri¬ 
bution  in  the  investigated  structure. 

3.2.1  Loading  and  Stress  Distribution 

In  general  the  forces  and  moments  acting  on  ,thc  column  of  a 
horizontal  boring  and  milling  machine  can  be  reduced  to  throe 
cases  of  loading  shown  in  Figure  3.26. 

a)  The  bending  moment  M  involves  bonding  deformations. 

b)  The  two  forces  F/2  acting  on  the  slide-ways  in  the  same 
direction  are  directly  introduced  to  the  sido  walls  of 
the  column.  They  cause  bending  deformations  duo  to  the 
moment  F-*  1  and  shear  deformations  duo  to  the  force 

Q  -  F. 

c)  The  torsional  moment  F  •  b  enn  be  reduced  to  the  couple 
of  forces,  which  is  applied  at  the  slide-ways  too.  The 
couple  of  forces  causes  twisting  of  the  cross-section 
and  under  certain  conditions  considerable  distortion 

of  the  cross-section  contour.  This  warping  leads  to 
considerable  twist  doformations  at  the  front  wall  of 
tho  column  and  by  this  at  tho  slido-wuys  for  the  cutter 
hond. 

Calculating  tho  olastio  behaviour  of  tho  specimen  following 
assumptions  concerning  tho  stress  distribution  are  mndo; 
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a)  The  wall  thle/.nos:;  is  wain  11  i n  relation  to  cross-suction 
and  length  of  the  structure.  Thor'eforu  normal  and  tangen¬ 
tial  stress  can  be  assumed  to  bo  constant  all  over  the  wall 
till  CkllUMM  , 

b)  Normal  stress  is  linear  along  the  walls  of  the  cross -sect ion. 
Tangential  stresses  are  distributed  ir.  the  form  of  a  ;)arnboln, 

c)  Along  the  walls  of  longitudinal  sections  tangent  ini  stresses 
are  assumed  to  be  constant, 

The  stiffness  of  the  walls  against  bonding  can  be  neglected  when 
regarding  the  direction  of  the  minimum  socond  moment  of  area.  In 
t ho  same  way  the  stiffness  of  the  walls  against  torsion  can  bo 
neglected. 

3. 2. 1.1  bonding  Doforma lions 

The  bending  deformation  of  a  hollow  bourn,  fixed  at  one  end  and 
loaded  as  shown  in  Figure  3.2G  a)  and  b),  can  bo  calculated  from 
the  well  known  formulas: 

a)  loaded  by  a  pure  bonding  moment  M 

f  —  ..  !~L P-  .-  /  2 

2EIy  L 

b)  loaded  by  a  single  force  y 

f _ f—t3 

'&  3  Ely  L 

3 . 2 . 1 . 2  Shear  bu format i ons 

’./lion  loading  a  beam  by  a  single  force  F,  not  only  normal  stresses 
but  also  shear  stresses  occur  in  the  cross-section.  They  lead 
to  additional  shear  deforma  I  ions,  The  deflection  In  the  height  1 


(3.37) 

(3.3b) 


L 


of  tlu*  1  oh  clod  n  roHM-Mnn  t  i  on  duo  to  shout*  force  cart  hu  wri  tten  n :  i 


Cl.  30) 


The  factor  takes  into  account  tho  distribution  of  shorn’  stress 
over  tlio  cross-section,  Fcippl  £c>J  gives  thu  followin';  equation 
for  the  factor  : 

(3.-10) 


*■  4-/T‘-dA 


For  the  shear  stress  it  can  be  written 

-  TW  JV  dA 

A 

and 


(3  .  -1 1  ) 


(3.-12) 


ly  u  t 

For  calculating  a  thin  walled  rectangular  croHS-Hoctian  the 
factor  (sue  f3 , 1,3  J  )  can  be  written  ns 


(3. -13) 


f 


In  this  equation  the  a  hour  stresses  of  tho  unloaded  walls  no,  2 
in  Figuro  3.20  b  huvu  boon  neglected,  as  the  wall-thickness  is 
small  in  comparison  to  .  tho  cross-section .  HoWuver  for  small  valuun 
of  the  fnctor  oC  "  h/b  tho  shear  stress  ill  tho  unloaded  walls  have 
to  bo  considered. 

Considering  the  loaded  walls  no,  2  in  Figuro  3.2(1  for  calculating 
the  influence  of  shour  force  tho  equation  to  compute  af  (son  fll  ,  KlJ  ) 
bocomos 

96  *  5^37^?  (2*3+  I0**'  IS*  +  5)  (3.4-1) 


I 


Equ,  (3,44)  shows  that,  tho  distribution  of  tho  shear  stresson 
in  tho  structure  depends  only  on  the  factors  i.o,  on  the  uhapo 
of  tho  croBM-Bootion.  Tho  relation-ship  botWuuiufand  both  tho 
factors  #’and  it  is  plottud  in  Figure  3,27.  The  diagram  uhowM  that 
tho  influonoe  of  the  shonr  stresses  in  the  unloaded  walla  can 
be  nogleated  for  or  ^  1  in  practical  cases  aa  tho  difforrnoe  bet¬ 
ween  at  and  at'  ia  small .  Only  in  the  range  of  <  1  these  tangential 
atroBBOB  hnvo  to  bo  considered. 

3,2.1.  3  Total  Deflect  1  on  duo  to  a  Single  Forco  F 

For  a  thin-Wallud  beam  loaded  by  a  Binglo  foroo  F  tho  total  do- 
floatioh  at  the  point  of  force  application  oan  bo  Written  an 

ftot  "  fU  +  fQ 

5  2 

With  the  rations  Q  ”  •  E  for  Htoel  and  C  -  ^  E  for  oast 

iron  and  using  tha  abbrovintion 

l/ll  -  14 


tho  total  deformation  of  tho  thin  wallud  buam  with  rectangular 
cross-suction  boooines  for 


a)  calculating  a  stool  structure: 

F( 3 


(3  +  <x)  W  *  (2 «  +  10*  *  +  (:) ’  45 } 

b)  ouloulating  a  oaut  iron  struoturu 

/>•«  '  ZiEth**)  </3>f  3 «•'*“) 


In  Figure  3,29  tho  deflection  duu  to  shear  force  is  compared  to 
the  deflection  due  to  bonding.  The  ratio  of  tho  dofloctions  at 
tho  point  of  forco  applico.tion  la  plotted  ovor  the  distanco  of 
tho  loaded  sootion  from  the  fixed  end  of  the  beam,  which  is  re¬ 
lated  to  tho  height  of  the  cross-section.  Tho  diagram  shows  that 
in  tho  range  of  usual  sizes  of  maohino  tool  columns 

-  L/h  -  2.5  +  4.0  with  1/L  >»  0. 15  +  0.0 

whore  L  moans  the  totnl  lenght  of  the  column  - 

tho  deflection  duo  to  tho  shear  force  has  considerable  signi¬ 
ficance  besides  tlio  bending  deflection.  Tho  influence  of  shear 
force  can  bo  neglected  in  comparison  to  bonding  deformation  only 
for  relatively  great  distances  of  tho  loadod  soction  from  the 
fixed  end  of  the  aolumn. 

It  may  be  assumed  for  instance  that  dofloctions  duo  to  shear 
force  can  be  negleated  if  their  part  of  the  whole  deflection 
is  less  than  15  %,  In  this  oase  the  influence  of  shear  remains 
disregarded  for  a  aolumn  with  a =  2  for  1  >  3  h.  For  a  beam 
with  Square  aross-section  shape  (oC-  1)  this  border  is  growing 
up  to  1  «4  h.  However,  if  the  column  has  a  small  height  of 
oross-Boatlon  in  comparison  to  its  breadth  b,  as  it  is  valid 
for  values  of  o c»  0,8  or  less,  shoar  deformation  can  only  be 
neglected  for  /?  >  G  (1  >  G  h). 

3 , 2 , 2 . 1  Torsional  Loading  and  Doformatlo n  of  Columns 

Tho  torsional  doformntion  of  a  column  fixed  nt  one  end  and 
loadod  as  shown  in  Figure  3.29  can  bo  calculated  by  well 
known  formulas  if  thoro  is  no 


distortion  of  tho  cross  section. 


Mr 

T17 


L 
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(3.47) 


For  thin  walled  columns  the  polar  area  moment  of  inertia 
usually  is  calculated  from  tho  equation  known  as  "Bredt's"- 
Formula 


A„  is  tho  area  enclosed  by  tho  aontro  lino  of  the  walla,  U 
tn 

is  the  oircumferonce  given  by  tho  center  linos  of  the  walla 
and  s  is  tho  variable  wall  thickness.  If  the  wall  thickness 
is  constant  around  the  circumference  the  oquation  for  calcu¬ 
lating  the  area  momont  of  inertia  becomes 

IT  -  (3.49) 

The  equations  (3.48)  and  (3.49)  suppose  tho  angles  remaining 
constant  under  loading  and  the  column  being  able  to  arohe  freely. 

Theoa equations  also  do  not  consider  arching  of  the  outside  walls 
and  distortion  of  the  cross-seotion.  The  torsional  loading 
has  to  be  applied  to  the  column  with  respect  to  the  flux  of 
shear  stress  being  aonstant  all  over  the  cross  soction  (T.j-  const,). 


As  tho  deflection  of  a  aolumn  caused  by  torsional  loading  due 
to  a  couple  of  forces  as  shown  in  Figure  3,29  cannot  be  compu¬ 
ted  by  normally  used  equations,  a  method  to  calculate  these 
torsional  deflections  has  been  evaluated  taking  into  account 
the  distribution  of  the  shear  stress.  Similar  to  bending  a  factor 
X  (see  f  3  ,  IS J  )  is  stated  which  considers  the  cross-sectional 
deformation.  Tho  index  of  x  corresponds  with  the  numbers  of  walls 
in  Figure  3,26 

-  -^T  (t2<X3+  2c(Zt  oC  +  0,2)  (3. SO) 


02  +  0,2 


(3.51) 
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3. 2. 2. 2 


Torsional  Deformations  with  and  without  Warping 
Effects 


For  the  pure  torsion  the  twisting  of  the  column  in  the  force 
noting  height  causes  a  deflection  of  the  walla  which  carry 
the  elide  ways.  This  deflection  can  be  calculated  from  al¬ 
ready  presented  equations 


SIr  ,l*n' 

The  torsional  moment  is  assumed  to  be 


,  UL 

b 


lip  -  F  .  b' 


and  with  the  area  moment  of  inertia 


I 


T 


4  Am  f 

D 


(3 . 4d  ) 


the  equation  to  aaloulato  tho  deformation  becomes 

f*  m  jS/m±^L 
’h  8  £s  <x 


(3.52) 


For  any  given  oroas-seotion  tho  relationship  botwoon  tho 
torsional  deformation  with  distortion  can  bo  writton  ns 


3<xs  f  2  a  f  0  5] 

s<i  +  <xFL  '/?/  J 


(3.53) 


In  tho  diagram  of  Figure  3. 30  tho  ratio  i’^  /f^*  is  plotted 
to  the  normalized  height  ft  -  1/h,  Tho  curvft  shfiws  that  in 
a  hoiGht  of  ft  -  2  the  torsional  deformation  due  to  a  distor¬ 
tion  of  the  croes-oootion  is  about  1.0  times  htchor  than  tho 
corresponding  torsional  deformation  without  distortion. 
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3.2.3 


Modol  Tests 


To  ohook  tho  investigation  methods  on  one  side  and  prove  the 
equations  derived  above  on  tho  other  side  experimental  model 
tests  have  been  carriod  out  using  Welded  Bteel  models. 

3. 2. 3.1  Tested  Columns  and  Teat  Method 

Tho  dimensions  of  the  wolded  Btool  models  have  boon  scheduled 
in  suoli  a  way  that  the  specimen  oan  be  handled  easyly  and  inve¬ 
stigated  With  os  small  expenditure  as  possible.  The  design  of 
tho  tested  columns  and  their  dimensions  aro  shown  in  Figure  3,31 
and  3.32  .  Tho  ratios  of  tho  dimensions  of  tho  modol  with  square 
orosB-aootioh 

o <*h/b  <*  1 

L/h  -  3.20 

s/h  ~  0,020 

corrospond  approximately  to  the  values  used  at  real  columns  of 
machine-tools. 

Tho  models  with  rectangular  cross-soction  have  tho  same  height 
L  -  10,7  In  (300  mm ) ,  as  tho  models  with  square  uross-seation 
have  but  tho  ratios 

Cf  -  2,00  or  cf  »  0.-10 

for  the  dimensions  of  tho  urosb-buotion.  All  models  aro  mnnu- 
faaturud  with  and  without  top  -  plato. 

Figuro  3,33  shows  the  toHt  rigs  to  investigate  tho  bohaviour 
of  tho  columns  under  bonding  and  torsional  loads.  Tho  columns 
aro  fixed  by  16  screws  on  to  a  stiff  intormodinte-plato  which 
itHolf  is  fautonud  by  screws  on  a  fixing-plato .  Tim  fixing  oon- 
ditlons  of  tho  columns  have  proved  very  stiff.  Novurtholuss 
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a  certain  flexibility  of  tho  clamping  causos  an  additional  da- 
flootion  of  the  column  and  thoroforo  tho  vortioal  deflections 
of  tho  flnngo  of  tho  column  aro  measured  to  separate  tho  de¬ 
flection  of  the  column, duo  to  tho  clumping  condition  from  the 
deformation  of  the  column. 

The  foraas  aro  appliod  to  the  column  by  moans  of  two  scrows  which 
arc  situated  in  tho  front  wall  aa  near  to  tho  aide  walla  as  possib¬ 
le.  For  bending  the  columns  aro  loaded  by  a  single  force  F  in  a 
way  that  two  forces  each  of  the  value  F/2  are  acting  on  the  side 
walls  in  tho  same  direction.  The  torsional  moment  represented  by 
a  aouple  of  forces  is  symmetrically  appliod  to  the  column  by 
moans  of  an  U-iron  that  is  loaded  by  two  hydraulic  load  units 
which  ara  fed  with  oil  of  the  same  pressure  by  ono  pump.  Apply¬ 
ing  tho  toraiondl  moment  it  is  ensured  that  ho  singlo  force 
occurs  which  could  lead  to  additional  bending  and  shear  deformation. 
Ill  three  different  heights  the  columns  are  loaded  by  bonding  foraos, 
or  torsional  moments,  tho  hoights  aroj 

oaso  at  height  of  loaded  section  1  «*  17.9  in  (405  nun) 

ottBo  bi  "  "  u  it  1  14.55  in  (370nun) 

oase  cl  1  -  0.3  in  (21.0  mm) 

All  deformations  and  deflections  aro  moasured  with  inductive 
deflection  pick-ups,  Tho  pick-ups  aro  mounted  on  a  stand  which 
is  fastened  on  tho  fixing  pinto  in  a  manner  that  the  stand 
does  not  move  whilst  loading  tho  column.  From  this  stand  tho 
abuoluto  deformations  can  bo  monsurod. 

3 . 2 . 3 . H  Investigation  of  Unribbod  Column ::i 

In  order  to  chock  tho  noouracy  ol’  tho  above  dorivod  methods 
to  calculate  tho  deformations  of  columns  under  bending  and 
torsional  loads  models  of  unribbod  columns  aro  tostod.  In  nckli- 
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tion  tho  bonding  stiffnoss  obtained  from  these  tosts  is  chosen 
as  a  standard  to  relate  the  measured  stiffnoss  of  ribbed  columns. 

3.2. 3. 3  Pending  Deflections 

The  table  in  Figure  3.3d  shows  tho  calculated  and  measured  do- 
floctiomof  tho  columns  with  square  cross-section.  The  measured 
values  given  in  tho  table  are  taken  as  the  average  values  of 
tho  deflections  which  have  been  measured  at  the  4  corners  of 
the  loaded  cross-section.  It  can  bo  seen  that  there  is  quite 
a  good  accordance  between  the  dalculated  and  measured  values 
of  deformation. 

The  deflections  of  the  columns  with  rectangular  cross-section 
measured  in  tho  same  way  are  compared  with  the  oaloulated cb- 
flections  in  Figure  3.3S.  The  column  with  top  plate  shows 
a  little  higher  stiffness.  The  errors  are  less  than  10  %, 

The  ooouring  errors  are  caused  by  the  simplifying  assumptions, 
which  have  been  made  for  the  calculation,  as  well  as  by  errors, 
mado  in  measurement,  during  running  tests,  and  manufacturing 
the  columns. 

3. 2. 3. 4  Torsional  Deflections 

The  measured  deformations  of  the  columns  under  the  applied 
torsional  moment  do  not  show  a  pure  torsional  deformation  but 
a  cross-sectional  distortion  which  is  drawn  in  Figure  3.38,  For 
this  reason  a  "stiffness  against  torsional  load"  is  deolared  as 
the  measured  deformation  at  the  point  of  loading  devided  by  the 
force  itself.  In  the  Figures  3.37  and  3,38  the  relationship  between 
the  measured  deformation  and  that  being  computed  neglecting  the 
warping  effects  are  given.  Figure  3.37  describes  tho  values  ob¬ 
tained  from  the  oolumns  with  square  cross-section  whilst  Figure  3,38 


shows  tHoso  of  the  column  w.i.  th  rectangular  cross-soc  L  ion .  Thu 
maasurod  results  n;:i’co  quite  well  with  tiiu  computed  ones  in 
tho  loading  heights  af  17.0  in  and  Id. 33.  The  accordance  becomes 
worse  in  n  height  of  0.3  in.  Those  differences  are  led  bn olt  on 
tho  influence  of  tho  flange  diminishing  the  cross-sectional 
distortion. 

3. 2. 3. 5  Investigation  of  nibbed  Columns 

In  order  to  compare  the  influence  of  some  kinds  of  internal 
ribbing  the  measured  deformations  of  tho  columns  loading  them 
in  a  height  of  17,9  in  are  only  taken  Into  account, 

3 . 2 , 3 , 0  Pending  stl Tineas 

Tho  inenaurod  bondi ng  stiff  nossoa  of  the  internally  ribbed  columns 
are  related  on  the  one  of  the  unribbod  column  ,  In  Figure  3,39 
the  normaliaod  bonding  stiffnesses  are  given.  The  listing  shows 
that  tho  stiffness  increases  by  ribbing  approximately  in  pro¬ 
portion  to  the  growth  of  the  area  moment  of  inertia.  Of  all 
tested  columns  tho  one  with  double  diagonal  ribbing  stiffens 
the  columns  best.  The  increase  of  stiffness  amounts  up  to  32  ',1 
whilst  tho  othor!  tested  typos  of  vortical  ribbing  increase  tho 
bending  stiffness  between  .1.3  %  to  21.  %  corresponding  to  the  kind 
of  ribbing,  Y/hen  considering  the  amount  of  material  in  the  ribs 
which  increases  the  bending  stiffness  by  relating  the  stiffness  to 
the  total  weight  of  each  column  re  poctivoly  it  becomes  obvious 
when  lo  ilcing  at  Figure  3.39  that  the  weight  increases  quicker  than 
the  stiffness.  Thus  tho  bonding  stiffness  of  tho  unribbed  column 
related  to  tho  weight  is  higher  than  those  of  the  internally 
ribbed  columns. 

3 . 2 , 3 , 7  ToraionaJ  Stiffness 

As  shown  in  Chapter  3. 2. 3. 4  tho  symmetrically  applied  couple  of 
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forces  distorts  tho  cross -section  of  tho  unribbad  columns 
£incl  thus  n  considerably  lower  torsional  stiffness  is  caused. 

For  this  reason  tho  calculated  stiffness  duo  to  a  pure  tor¬ 
sional  deformation  of  tho  column  has  boon  chosen  ns  a  base  to 
relate  tlio  measured  stiffnesses  of  tho  ribbed  column.  As  it  is 
a  theoretically  calculated  stiff ness  this  valuo  is  marked  in 
Figure  3.15  as  "ttioorctical" .  For  comparison  only  tho  tost  re¬ 
sults  of  the  loading  height  of  17,9  in  aio  tnkon  into  adoount. 

Tho  normalized  stiffnesses  against  torsion  are  listed  in  Figure 
3 .  -10.  in  opposite  to  tho  bonding  stiffness  the  stiffness  against 
torsion  can  bo  increased  by  a  top-plato  which  becomes  ovidont 
comparing  the  unfibbod  and  the  ribbed  columns  with  parallel  ribbing. 
Tho  stiffness  against  torsion  of  those  columns  without  top-plate 
varios  between  0.G  %  and  15  %  at  the  maximum  whilst  tho  columns 
with  top-plata  reach  51.5  %  of  tho  theoretical  stiffness.  With  ono 
diagonal  rib  tho  column  reaches  a  stiffness  of  05  %  and  tho 
highest  stiffness  is  measured  at  the  column  With  double  diagonal 
ribbing,  and  amounts  up  to  1.15  of  the  theoretical  stiffness.  The 
diagonally  ribbed  columns  with  top-plato  even  shew  a  higher  stiff¬ 
ness  but  tho  difforonao  to  tho  ones  without  top  plate  is  not  as 
important  as  at  tho  unribbed  oolumns  or  those  with  ribs  parallel 
to  the  outor  walla.  Also  good  results  show  tho  columns  with 
horizontal  ribs,  those  ribs  suppress  the  cross-moot ional  distor¬ 
tion  and  thus  stiffen  the  columns  up  to  a  stiffness  equal  to 
the  stiffness  of  pure  torsion. 

The  torsional  stiffnesses  related  to  tho  total  weight  of  tho  columns 
show  a  tendency  similar  to  the  one  of  tho  bonding  stiffnesses, 
Choosing  the  torsional  stiffness  without  cross-sectional  distor¬ 
tion  as  a  baso  this  stiffness  is  not  reached  by  any  of  tho  inter¬ 
nally  ribbed  columns,  From  Figure  3.40  can  bo  seen  that  tho  effec¬ 
tiveness  of  the  material  in  ribs  is  best  whoa  ribbing  tho  columns 
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with  horiaontnl  fibs,  two  horizontal  ribs  and  n  top-plnto  qivo 
«  torslionnl  utiffnosei  minted  to  tho  v/oi tsh L  of  OB  %  of  the  theo¬ 
retical  stiffness, 

3,3.4  Teat  Method  of  Dynamic  Invostlpatlonu 

In  order  to  obtain  characteristic  data  do  tonal  ninij  tho  dynntnlo 
bohaviour  the  columns  iiavu  boon  dynamically  tested,  They  nro 
loaded  by  a  dynamic!  force  of  +  100  p  which  in  (juliornted  by  an 
oleotro-dy numic  oxcitor,  Tho  experimental  urranttomuht  la  shown 
ill  tho  blaokdiuf;rnm  of  Figure  3,41.  According  to  bonding  and 
torsional  loads  tho  dynamic  force  is  applied  at  three*  different 
points  to  tho  column  us  shown  In  Figaro  D ,  43,  tho  pick-ups  aro 
mounted  on  the  column  opposite  to  the  oxditcr  to  measure  tho 
response  aurves,  iiixoiting  tho  columns  at  a  Corner  on  tho  top  tho 
columns  vibrato  at  different  frequencies  in  their  bunding  and 
torsional  modes,  a  few  plate  vibrations  are  incited  as  Well, 
Applylnu  tho  exciter  force  in  tho  centre-line  on  thu  tup  of 
t’lu?  ool.UMrtfl  tiie  bonclMuj  modo  and  pluto  vibrations  nru  inoltud 
and  when  oxoltLuii  tho  columns  in  tho  uontro  of  a  sldo  wall 
mainly  pinto  vibrations  occur, 

In  all  throe  positions  of  the-  uxoi tor  tho  response  curves 
lir.vu  picked  up  inn  frequency  range  botwuon  UOO  ops  and  3000  ops, 
Thu  body  moden  and  pinto  vibrations  am  identified  by  mouHUrlng 
tho  amplitudes  and  dirooliouu  of  the  vibrations  at  many  points  of 
tho column  so  that  the  modus  can  bo  plotted.  Of  the  body  modes  and 
and  of  iiomo  plate  vibrations  the*  dampiin:  coefficients  have  boon 
detorminnd  from  thu  decay  curves, 


3.3.4, 3 


ltoapoiwy  Curves  and  Mode  beg  pcs  of  Unrlbbod  and 
i nturnally  nibbed  Columns 


111  Figure  3,43  tlio  renonance  curves  of  U  columns  without  top- 
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pinto  wlion  exciting  tho  columnu  at  tv  corner  on  tho  top  nro  uhown. 
Tho  roaonanou  curves  differ  in  thoir  shape  oxoopt  that  two  re¬ 
sonance  punks  can  bo  found  at  almost  tho  name  frequencies  •  Tlioy 
m’o  the  rosonnnca  l’raquonciou  of  tho  bonding  modo  nt  450  opa  and 
tho  ono  of  tho  torsional  modo  at  1300  opa.  These  resonance  fre¬ 
quencies  change  only  vory  littlo  although  tho  stiffnessou  of  tho 
columns  ngninut  bending  and  torsional  loads  altoi  for  relatively 
largo  nmottnt  corresponding  to  tho  difforont  typQB  of  ribbing.  But 
us  tho  vibrating  mnutt  id  changed  by  tho  riba  aa  Woll  dtiffnoaB 
and  mass  arc  altered  about  the  uunie  mount  and  tltUH  tho  resonances 
of  tlio  bending  nnd  torsional  modoa  remain  at  thoir  i’roquenoioB 
ihdopahdont  from  tho  typo  of  ribbing, 

A  resonance  frequency  Which  Changes  between  450  ops  and  750  Opa 
according  to  tho  difforont  internal  ribbing  iu  tho  ono  of  tho 
toraiotuvl  modo  with  orohis-ucatiuttnl  deformation.  The  frequencies 
of  thin  modo  shape  vary  in  such  a  largo  range  because  the  stiffness 
of  tho  oroHS-Hoction  against  distortion  alters  ho  ltiuah, 

The  modo  shapes  of  those  M  body  modes  of  the  whole  column  are 
drawn  in  Figure  3.44.  They  are  tho  modo  shapoS  Which  are  of  in¬ 
tercut  in  a  dynamic  analysis  of  the  ntruaturo. 

Any  resonances  else  which  have  boun  measured  at  tho  columns 
have  boon  dotormil, jd  as  plate  vibrations  of  tho  outer  walls  or 
of  the  ribs  respectively.  [Jut  unlike  to  tho  resonances  of  tho 
body  modus  no  relationship  can  bo  found  between  tho  rosonanoe 
froquohuios  and  mode  shapes  of  tho  pinto  vibrations'  in  the 
roHponuo  curves  of  tho  difforont  columns.  Tho  plate  vlbrutioiui 
are  ini’ltioticod  by  the  internal  ribbing  and  thus  tho  ruHonnnoo 
ire  quo  no  ion  of  uimilur  shape  of  tho  pinto  vibrations  differ  but- 
woon  tho  columnu. 

Most  resonance  fruquoncius  of  pluto  vibrations  nave  bouiv  deter- 
ininud  at  the  unribbud  uud  diagonally  ribbed  columnu.  dome  of 
those  modo  shapes  arc  represented  in  I'M gvru  h.4b.  At  tho  columnu 
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with  horizontal  riba  and  those  with  internal  riba  pantile!  to 
the  outer  walla  the  number  of  resonances  of  plate  vibrations  is 
smell  in  the  froquenoy  range  up  to  300Q  eps.  Thin  is  led  back  to 
the  fact  that  the  horizontal  riba  and  the  or. os  prallel  to  the  out¬ 
er  walla  do vi do  them  into  smaller  diaphragms.  Thus  resonances  of 
lowor  frequencies  are  suppressed, 

Not  only  internal  ribbing  but  also  a  top-plate  which  closes  the 
column  at  the  top  diminishes  tho  number  of  plate  vibrations  as 
it  can  be  seen  from  a  comparison  of  the  response  curves  of  the 
columns  with  and  Without  top-plate.  These  response  curves  are  re¬ 
presented  in  Figure  3,40.  At  the  aolumn  With  top-plate  even  tha 
torsional  mode  With  cross-sectional  distortion  does  not  appear, 
because  the  exciting  force  is  in  thio  case  applied  to  the  aolumn 
near  to  the  top-plate  and  tile  top-plate  suppresses  tin  eroraS-rjoo- 
tional  distortion i 

For  two  other  points  of  force  application  the  ronpomio  curves 
of  tho  0  columns  aro  shown  in  Figure  3,-17  and  3,40.  The  results 
derived  from  thedo  curves  are  tho  ones  that  for  improving  tho 
dynamic  behaviour  of  oolUmnB  plate-vibrations  must  bo  suppressed 
as  it  iu  possible  for  example  With  horizontal  ribs.  Another  way 
to  avoid  plate  vibrations  is  ta  stiffen  the  plates  by  wall  rib¬ 
bing.  In  order  to  oheolc  the  influence  of  wall  ribbings  on  tho  dy¬ 
namic  behaviour  of  aolumus  somo  typora  of  wall  ribbing  huVo  boon 
tested,  'file  results  of  theae  tost  are  discussed  in  the  following 
ohaptoro. 


3 , 2 , U  Touted  Columns  with  Wall  lUbbing 

Tho  dimensions  of  tho  columns  with  wall  ribbing  arc  like 
those  of  the  internally  ribbed  columns  and  they  are  touted 
in  an  equal  way.  In  this  investigation  throe  different  types 
of  wall  ribbing  are  tested,  They  are  shown  in  Figure  3,49. 
The  columns  of  type  A,D  havo  horizontal  wall  ribs  around 
tho  column,  Typo  U,K  iu  derived  from  typo  A  by  adding  a 
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vortical  i'ib  ill  tho  centre  lino  of  each  wall  from  the 
bol Low  to  tho  top  of  tho  column  and  typo  0,1'’  ia  a  variation 
of  typo  A  which  is  originated  be  dividinG  the  fields  between 
tho  horizontal  ribs  diagonally, 

3,2, 0,1  Itegults  of  Static  Testa 

An  tho  Wall  ribs  add  only  a  little  to  tho  area  momont  of  inertia 
but  on  tho  other  hand  tho  weldings  Weaken  the  structure  be¬ 
cause  of  the  notch  effect  the  columns  With  wall  ribbing  do  not 
roach  tho  bending  stiffness  of  tho  unribbed  column  which  bo- 
cometi  obvious  by  looking  at  tho  diagram  in  Figure  3,50, 

Tho  stiffness  ugainat  torsion  yot  has  been  improved  a  little 
by  the  wall  ribs.  This  mainly  is  cauaod  by  the  horizontal 
ribs  which  detain  the  orujs-aectional  distortion.  The  tor¬ 
sional  stiffness  of  columns  with  wall  ribs  is  listed  in 
Figure  3.01.  From  that  figure  can  be  seen  that  the  columns 
of  type  F  with  4  horizontal  and  diagonal  wall  ribs  is  stiffost 
although  there  still  remains  a  remarkable  amount  of  cross- 
soational  distortion. 

3 , 2 , G , 2  Results  of  Dynamic  Tests 

Figure  3,02  and  3,03.  show  response  curves  of  the  tested  columns. 
In  Figuro  3.82  response  aurves  of  six  columns  are  given  when 
exulting  them  at  a  corner  on  the  top.  Those  response  curves 
indioato  that  the  wall  ribs  do  not  change  tho  frequencies  of 
the  body  inodes  but  avoid  a  lot  of  the  plate  vibrations  which 
aro  incited  at  tho  unribbod  column  (soe  Figuro  140 .  Looking  at 
Figuro  3,33,  which  shows  the  responue  curves  of  the  columns 
when  exciting  them  on  the  top  in  the  centre  of  a  wall  the 
el feotivenesa  of  the  wall  ribs  booomos  obvious.  Up  to  frequen¬ 
cies  near  to  the  one  of  the  pure  torsional  mode  pinto  vibrations 
can  bu  avoided  by  wall  ribbing,  In  dynamic  tests  of  large  machine 
tools  it  has  been  proved  that  the  frequency  of  tho  pure  torsional 


mods  of  the  column  can  be  assumed  as  a  border  up  to  which  reson¬ 
ance  frequencies  of  columns  are  of  Interest,  Frequencies  of  n 
column  higher  than  the  one  of  the  torsional  mode  In  most  oases 
do  not  influence  the  dynamic  characteristic  of  the  whole  machine 

3,2,6  Damping  Coefficients  and  Mode  Shapes 

Besides  the  number  of  resonances  and  the  resonance  frequencies 
the  damping  coefficient  is  an  Important  parameter  whloh  deter¬ 
mines  the  dynamia  characteristic  of  a  structure.  As  it  is 
necessary  to  know  the  exact  damping  coefficients  for  a  calcu¬ 
lation  of  cross  response  curves  the  damping  coefficients  in 
the  different  mode  shapes  of  the  tested  columns  have  been 
measured.  The  results  of  the  damping  measurement  are  plotted 
in  the  diagram  of  Figure  3,64,  From  the  diagram  can  be  seen, 
that  for  different  body  mode  shapes  different  damping  coeffi¬ 
cients  exist.  In  the  bending  mode  the  highest  damping  coeffi¬ 
cients  have  been  measured.  The  coefficients  vary  in  a  large 
range  thus  no  definite  value  can  be  given  for  a  computation, 

For  the  torsional  mode  with  cross-scotional  distortion  the 
damping  coefficients  lie  near  a  mean  value  which  amounts 
a  quarter  of  the  mean  value  of  the  one  of  the  bending  mode, 

The  damping  coefficient  of  the  pure  torsional  mode  is  a  little 
higher  and  varies  a  bit  wider  than  the  ones  of  the  torsional 
modo  with  cross-sectional  distortion. 

From  Rg.  3.  54  can  be  deduced  that  for  each  body  mode  a  damping 
coefficient  exists  in  a  certain  range  which  differs  definitely 
from  the  others.  Although  at  the  moment  it  is  not  possible 
to  prodict  the  damping  coefficients  of  a  structure  it  Is 
thought  possible  that  by  systematic  Investigations  a  list 
of  damping  coefficients  can  be  set  up  in  order  to  enable 
more  exact  computations  of  cross-response  curves. 

The  plate  vibrations  in  most  cases  are  not  of  interest  in 
a  dynamic  nnalysis  of  a  machine  tool  structure  because  in 
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ft  well  designed  structure  plate  vibrations  must  not  occur. 

Tims  only  ft  few  damping  coefficients  of  different  modes  of 
plato  vibrations  hnvo  boon  measured  to  chock  the  amount  of 
those  coefficients .  They  are  all  near  the  damping  coeffi¬ 
cients  of  the  torsional  modes. 

3,2.7  Deductions  from  the  Modol  Tests 

From  the  static  and  dynamic  invos tigntions  of  the  simplified 
models  of  columns  a  few  important  points  of  view  can  be  de¬ 
rived  which  should  bo  taken  into  account  when  designing 
columns  of  machine  tools. 

To  obtain  a  rospeotablo  statia  stiffness  of  a  column  the  stiff¬ 
ness  against  bending  and  torsional  loads  have  to  be  determined 
independently  because  a  structure  whloh  is  stiff  against  ben¬ 
ding  is  not  necosr.arily  stiff  against  torsional  loads.  To 
optimizo  the  bonding  stiffnosa  of  a  column  the  cross-seotion 
should  be  designed  ns  large  as  possible  or  for  given  dimen¬ 
sions  of  a  cross-section  the  outer  Walls  should  bo  thickened 
instead  of  adding  internal  ribs  to  stiffen  the  struoturo.  Any 
internal  ribbing  will  not  increase  the  bonding  stiffness  as 
much  as  an  oqual  amount  of  material  whioh  is  used  to  thicken  the 
outer  walla. 

In  Figure  3,55  a  list  of  normalized  cross-sectional  values  as 
area,  aroa  moment  of  inertia  and  area  moment  of  inertia  relatoci 
to  tho  area  of  unribbod  and  ribbed  cross-sections  is  given. 

Ail  cross-sectional  valuos  are  related  to  tho  ones  of  an  un- 
ribbed  square  cross-section  6. 2d  in  by  6.24  in  and  0,157  in 
(4mm)  wall  thickness,  For  each  ribbod  cross-section  the  afore 
mentioned  valuos  of  two  unribbod  cross-sections  are  computed 
when  distributing  tho  material  of  tho  ribs  evenly  to  tho  outer 
walls.  In  the  first,  caso  tho  wall  thickness  is  growing  to  tho 
innor  and  in  the  second  case  it  is  growing  to  the  outer  of 
tho  eross-soctJ on.  From  Figure  3.55  can  be  seen  that  i.e. 
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butins  the  material  of  tho  ribs  to  the*  outer  walls  the  area 
moment  (jrowu  by  09  %  or  74  %  respectively, 

Thu  torsional  stiffness  when  loading  u  column  by  a  couplo  of 
foroos  oaa  essentially  bo  inuroasod  by  internal  riba  which 
suppress  a  crosn-aoctional  distortion,  Anothor  Way  to  obtain 
a  aUfficiont  torsional  stiffness  id  to  choose  a  crosB-SoCtion 
Which  does  not  distort  if  tho  aotina  forcos  nro  applied  pro¬ 
perly.  Such  a  urosd-soation  for  example;  is  a  oirolo. 


Co  1 10  lUa  lulls 

With  regard  to  optimum  use  of  machine  tools  in  manufacturing 

plants  it  is  important  ntul  necessary  to  know  those  machining 
conditions  definitely!  where  the  installed  machine  power  can 
he  turned  to  profit  us  far  ns  possible  without  chatter  ocouring. 

After  having  extensively  studied  the  fundamentals  of  the  chat¬ 
ter  behavior  of  machine  tools,  n  procedure  could  be  developed, 
and  has  been  detailed  in  thiH  report,  by  menus  of  which  clear 
and  understandable  diagrams  of  tin*  chatter  behavior  of  latheH 
nml  milling  machines  can  be  produced  theoretically,  based  upon 
the  theory  of  regenerative  chatter  a  mathematical  model  of  the 
system  machine  tool  -  cutting  proaesB  has  been  established,  of 
which  the  stability  analysis  can  easily  bo  onrried  out  by  apply¬ 
ing  the  NyquiHt-Criteriou, 

In  prnotioe  numerous  parameters  and  mnohining  conditions  must  be 
taken  into  aoooUnt t  when  carrying  out  a  complete  stability  ana¬ 
lysis  of  a  machihe.  Because  of  this  fact  different  digital 
computer  programs  have  been  developed  for  the  quick  and  economi¬ 
cal  performance  of  extensive  chatter  investigations  in  consider¬ 
ation  of  t lie  practically  interesting  conditions.  Some  results  of 
shatter  investigations  on  a  lathe  and  a  milling  machine  have 
been  disoUHsud  ns  oharao ter 1st io  examples, 

From  theoretical  and  practical  studies  into  the  chatter  beha¬ 
vior  carried  out  so  far,  three  test  procedures  result,  which 
in  tho  future  can  oiler  the  possibility  to  provide  the  necessary 
data  for  specifying  maohlnu  tool  chatter  behavior.  However, 
further  extensive  investigations  and  discussions  between  the 
manufacturers  and  tho  muchino  tool  users  seem  to  bo  necossary, 
before  one  of  these  procedures  will  prove  ns  the  optimum  one 
and  before  definite  machine  tool  chatter  sped  Non  I.  ions  can  bo 
laid  down, 


The  methods  for  proenlculat ing  the  static  and  dynamic  behavior 
of  machine  tool  systems  «s  hriofly  done ri bod  in  this  linal 
Report  and  in  more  detail  in  the  according  Qviartorly  Reports 
are  found  out  vory  helpful  during  the  design  processes. 

These  methods  developed  under  this  contrnot  cover  n  wide 
range  of  problems  which  could  not  or  only  unsufficiently  ho 
solved  so  fnr. 

However,  it  should  be  pointed  ndt  that  still  problems  occur, 
when  dealing  with  the  statio  and  dynnmio  behavior  of  connection 
points  within  a  struaturee.  g,  between  elements  or  at  the  foun¬ 
dation  . 

If  these  parameters  -  as  In  several  cases  -  may  be  neglected, 
the  above  mentioned  methods  and  corresponding  programs  work 
very  Well,  In  oases  of  high  influences  duo  to  connection  - 
conditions  between  the  elements  of  a  system  the  overall  be¬ 
havior  oah  only  be  determinated  qUalitativly  so  far.  in  this 
oases  nevertheless  the  possibility  to  oompnre  some  given  sug¬ 
gestions  in  form  of  different  designed  structures  is  given  by 
dominating  the  above  mentioned  influences. 

In  a  thooretioal  analysis  some  equations  have  boon  derived 
to  calculate  the  bending  deflections  of  a  column  when  tuking 
the  influence  of  shear  force  into  account  and  moreover  equa¬ 
tions  have  been  formulated  to  caloulatu  the  torsional  defor¬ 
mation  of  a  column  loaded  by  a  couple  of  forces  which  distort 
the  croaa-Hootion .  The  influence  of  shear  force  and  the  one 
of  cross-sectional  distortion  ns  a  tunc lion  of  the  loading 
has  been  pointed  out . 

In  testa  of  models  of  uiiribbcd  columns  the  theoretically  de¬ 
rived  equations  havo  been  proved.  In  order  to  analyse  the  in¬ 
fluence  of  intornul  ribs  on  the  static  stlffnoss  a  number  of 
models  of  columns  havo  boon  tostod  and  tho  results  presented 
show  that  for  optimizing  the  bending  stiffness  the  material  is 
best  used  not  by  adding  ribs  mo  n  column  but  in  thickening 
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tbo  outer  walla  or  In  enlarging  the  cross-section.  Whilst 
the  torsional  stiffness  loading  a  column  of  rectangular  cross- 
section  by  a  couple  of  forces  is  increased  in  the  best  way 
by  adding  diagonal  or  horizontal  ribs  or  by  choosing  a  oircular 
oross-seot ion .  Furthermore  it  turned  out  from  the  statio  tests 
of  columns  with  wall  ribs  that  this  kind  of  ribbing  does  not 
sufficiently  stiffen  the  aolumns  neither  against  bending  nor 
ugainst  torsional  loads. 

Plato  vibrations  can  create  relatively  largo  deformations  of 
n  column.  From  dynamic  tests  of  the  internally  ribbed  columns 
it  turned  out  that  plate  vibrations  can  be  avoided  at  lower 
frequencies  by  dividing  the  large  areas  of  the  side  walls  by 
for  example  horizontal  ribs.  But  the  best  way  to  suppress  plate 
vibrations  is  thought  in  stiffening  the  walls  by  typos  of  wall 
ribbing  which  divide  the  walls  into  sufficiently  small  areas 
between  the  ribs.  This  postulation  has  been  proved  by  tests  of 
columns  with  different  types  of  wall  ribs.  The  dynamia  test  of 
these  columns  indicate  the  plate  vibrations  can  be  avoided 
by  wall  ribbing  up  to  frequonoies  near  to  the  first  pure  tor¬ 
sional  mode. 
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Figure  2.1  :  Influences  upon  the  Chatter  Behavior  of  Machine  Tools. 
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Fig,  3,22  Exciting  Scheme 
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